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THEORETICAL STUDY OE CORRUGATED PLATES: 


SHEARING OF A TRAPEZOIDALLY CORRUGATED PLATE WITH 
TROUGH LINES PERMITTED TO CURVE 


By Chuan-jui Lin* and Charles Libove** 
Syracuse University 


SUMMARY 

A theoretical analysis is presented of the elastic shearing of a 
trapezoidally corrugated plate with discrete attachments at the ends 
of the corrugations. Numerical results on effective shear stiffness, 
stresses, and displacements are presented for selected geometries and 
end-attachment conditions. It is shown that the frame-like deformation 
of the cross-sections, which results from the absence of continuous end 
attachments, can lead to large transverse bendir® stresses and large 
reductions in shearing stiffness. 


INTRODUCTION 

In a previous report (ref. 1) a theoretical analysis was presented 
of the elastic shearing of a trapezoidally corrugated plate with discrete 
attachment at the corrugation ends on the assumption that the trough lines (mn 
in fig. 1) are held straight. This assumption limited the applicability 
of reference 1 mainly to the case in which the corrugated plate is 
attached to a flat plate along its trough lines. 

The present report analyzes again the elastic shearing of a 
trapezoidally corrugated plate but removes the assumption that the trough 
lines are held straight. Thus the present analysis is applicable to a 
corrugated plate alone, rather than to a corrugated plate which is 
attached to a flat plate. The removal of this constraint reduces, of 
course, the effective shearing stiffness and alters the nature of the 
stresses and displacements. 

The analysis of reference 1 considered two kinds of conditions along 
the trough lines: (a) complete freedom of rotation and (b) complete 

suppression of rotation. In the present analysis only the first of these 
conditions is considered, as that is the only meaningful condition for a 
corrugated plate alone. 


*NDEA Fellow 


**Professor of Mechanical and Aerospace Engineering 


As in reference 1, three different types of discrete attachment at 
the ends of the corrugation are considered. These are illustrated in 
figure 2 and may be described as follows: 

(i) Point attachments at the ends of the trough lines only (fig. 

2(a)), the attachments being considered as mathematical points^ 
providing testraint against displacement but not against 
rotation. 

(ii) Point attachments at the ends of. both the trough lines and 
thd crest lines (fig. 2(b)), the attachments again being 
considered as mathematical points. 

(iii) Very wide attachments at the ends of the trough lines only, as 
Shown in figure 2(c). This kind of attachment is approximated 
in the analysis by means of the idealization shown in figure 2(d), 
i.e. by adding, to the end constraints of figure 2(a), end 
constraints against vertical displacement (but not against 
longitudinal displacement) at the junctions of the trough plate 
elements and the adjacent sloping plate elements. 

In cases (i) and (ii) no consideration is given to the possibility that the 
member to which a corrugation end is attached will interfere with the 
deformation of the corrugation. 


The analysis is based on the method of stationary total potential 
energy. Each cross section is assumed to have certain degrees of freedom 
for deformation in and out of the plane of the cross section. By equating 
to zero the first variation of the total potential energy, differential 
equations and boundary conditions are obtained for these degrees of freedom 
as functions of the longitudinal coordinate (z in fig. 1). Solution of 
these equations leads to all the desired information. 

Numerical results on shearing stiffness, stresses and deformations 
for selected families of geometries are presented and discussed. 
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SYMBOLS 


A ll* A 12’ etC * 


A l* A 3 


A !» A 4» A s» A 8 


A*, A*. A*. A* 

A l’ A 4’ A 5’ A 8 
A 1 * A 2 ’ A 3 * A 4 


A 6 


A ll’ A 22 ’ A 12 


a l* a 2 * a 3 


a ll* a 12’ a 22 


a *l’ a *2 ’ a 22 


a il’ a l2 ’ a 22 


22 


a* 

22 


22 


a 22 


Bj/A (j=l,2,..,8) 


B*. B* 
fi- 


de fined by equations (13b) 

coefficients in the displacement equations (C39) 

coefficients in equations (B50a) through (B50e) 
for displacements; obtained by solving equations 
(B52) , (B53) or (B54) , depending on the type of 
attachments at the ends of the corrugations 

coefficients in equations (B42) for displacements 

defined by the first four of equations (30) 

coefficients in equations (C33) for displacements 

coefficient in displacement equation (C39) 

defined by equations (C8b) 

defined by equations (40) 

coefficients in expression for U, (see eqs. 

(11) and (12)) 

defined by equation (B7) 

defined by equations (B25) 

defined by equations (C7) 

defined by the first of equations (C14) 

defined by the first of equations (C22) 

defined by equation (D6) 

defined by the first of equations (D17) 

obtained by solving equations (B23) 

coefficients in equations (D34) for displacements 

coefficient in equations (D34) for displacements 
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b 



b Jl’ b 12’ b 22 


C j /A j (j=l,2, . . . ,8) 


c oo* c ox’ etc- 


c* . . c* 

11 * 12 ’ 22 


C U’ C X 

r\ r\ 
cv, cy 


'00 


Gt/k 


D 


Dj/A^ (j=l,2, ... ,8) 

Dj/B (j=l»2,3,4) 
Dj/Aj (j=l,2,3,4) 

d oo* d ii» etc - 


d 10* d ll» d 21» etC - 


one-half th'e length of the corrugations 
(see fig. 1(b)) 

defined by equations (40) 

coefficients in expression for U t (see 
eqs. (2) and (3)) 

defined by equation (B7) 

defined by equations (B26) 

obtained by solving equations (B23) 

characteristic length (taken as pitch p in 
numerical work) 

coefficients in expression for U , (see eqs. 
(5) and (6a)) 

defined by equation (B7) 
defined by equations (B27) 
defined by equations (32) 
defined by equations (B60) 

frame flexural stiffness; see equation (14) 

obtained by solving equations (B23) 

defined by equations (D23) and (D24) 

defined by equations (C28) and (C29) 

coefficients in expression for Ug^ (see eqs. 
(5) and (6b)) 

defined by equations (B28) 
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d ll* d 21 


d 1> d 2 , etc. 


defined by equations (B62) 
defined by equations (B8) 

Young's modulus associated with frame bending 
of the cross sections 


Young's modulus associated with longitudinal 
extension 


J 10 


coefficient in equations (B42) , (B50a) , (B50b) 
and (B50e) for displacements 


10 


defined by the last of equations (30) 

one-half the width of the trough plate element 
(see fig. 1(a)) 


e oo* 

e ll’ 

etc. 

defined by equations (A2) 


e ll’ 

*12’ 

e 22 

coefficients in expression for U tw 
(9) and (10)) 

(see eqs. 

e oo’ 

*11’ 

etc. 

coefficient in expression for U s jj 
(5) and (6c)) 

(see eqs. 

e 



defined by equation (B7) 


e oo 



defined by equation (D5) 


*ii’ 

e *2 ’ 

e* 

22 

defined by equations (B30) 


e** . 
00’ 

e** , 
11* 

etc. 

defined by equations (B29) 


e** 

00 



defined by equation (D17) 


e ir 

e l2 ’ 

®22 * e l0’ e 20 

defined by equations (B62) 


p' ' 

e ll’ 

e i2’ 

e' ' 
22 

defined by equations (C6) 


0 * 
00’ 

; io* 

P* P* 0 * 

11’ 12’ 22 

defined by equations (CIO) 


e 22 



defined by equation (C14) 
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e 00’ *10’ e 22 


F 

f 

f 


f , 

f , 

etc. 

ss 

sc 


f , 

f , 

etc. 

ss 

sc* 



G' 

g 2i» 8 22 » ® 3 i ’ S 32 
g 
h 
h 

1 = vCT 


^1* J 2’ ^3 


3 

k 


k Q , k 2> etc. 


k 02* k 20’ etC * 


k 02’ k 20 * etC ’ 


defined by equations (C22) 
shear force (see fig. 1(b)) 

width of the crest plate element (see fig. 1(a)) 

defined by equation (B7) 

functions of z defined by equations (31) 

functions of z defined by equations (31) 

shear modulus associated with middle surface 
Bhear of the plate elements 

shear modulus associated with torsion of the 
plate elements 

defined by equations (39) 

defined by equation (B7) 

height of corrugation (see fig. 1(a)) 

defined by equation (B7) 

torsion constants of plate elements 01, 12, 23 
respectively (see eqs. (7) and (8)) 

defined by equation (B7) 

width of the inclined plate element (see fig. 1(a)) 

coefficients in characteristic equations (B12) 
and (B20) ; defined by equations (B13) 

defined by equation (B7) 

defined by equations (C21) 

defined by equations (D16) 
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L ll* L 12* etC * 


defined by equations (B24) 

^iji ) ^22 * etc - 


defined by equations (C27) 

L ll* L 12’ &tC ‘ 


defined by equations (D22) 

m 


defined by equation (B7) 

*4* N 5 


defined by equations (B56) 

*11* N 12’ etC - 


defined by equations (B55a) through (B55d) 

? 5 


defined by equation (B58) 




*41* *42’ *43’ 

*44 

defined by equations (B57a) 




*51 * *52’ *53’ 

*54* *55 

defined by equations (B57b) 

*41 5 *42’ *43’ 

*44 

defined by equations (B59) 

» ^2 5 


defined by equations (B9) 

p B , p c , p d , p e 


real numbers defined by equations (B46) through 
(B49); obtained by solving equations (B23) with 
j = 1 and 5 and noting equations (B31) 

P 3 


defined by equation (C41b) 

P ll* ? 12 ’ etC - 


defined by equations (C41a) 

P 


pitch of corrugation (see fig. 1(a)) 

p’ 


developed width of one corrugation, 2e + 2k + f 

Pi’ P 2 * P 8 

defined by equations (BIO) 

q b , q c , q d , q e 


real numbers defined by equations (B46) through 
(B49); obtained by solving equations (B23) with 
j = 1 and 5 and noting equations (B31) 

Q 3 


defined by equation (D37) 

Qj 2 * e ^c. 


defined by equations (D36) 
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q l* q 2* q 3’ q 4 
R 

R j (j-l,2,...,10) 

R ! 

R* 

R 

Rj (j*l,2 6) 

R 

Rj (j**l»2 6) 

r - R/c 

Tj H R^/C (j=l,2 
r «= R/c 


defined by equations (Bll) 

variable in characteristic equation (B20) 

roots of characteristic equation (B20) 

complex conjugate of R^ 

complex conjugate of R^ 

variable in characteristic equation (C19) 
roots of characteristic equation (C19) 
variable in characteristic equation (D14) 
roots of characteristic equation (D14) 
variable in characteristic equation (B12) 


r = R/c 


S 


B 


> 




real numbers defined by equations (B46) through 
(B49); obtained by solving equations (B23) with 
j = 1 and 5 and noting equations (B31) 


s r s 2 . 


transverse coordinates along the cross-sectional 
centerline (see fig. 3(a)) 


r\ 


sv. 


sy 


defined by equations (B60) 


TPE 

t b , t c , T°, t e 


r\ 

tu. 


r\ 

tx 


total potential energy of a single corrugation 


real numbers defined by equations (B46) through 
(B49) ; obtained by solving equations (B23) with 
j * 1 and 5 and noting equations (B31) 

thickness of corrugation (see fig. 1(a)) 

defined by equations (B60) 
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strain energy of an entire corrugation (see 
eq. (15)); also real part of and R 2 when 
R^ and R 2 are complex (see eqs. (B44)) 


strain energy per unit length of corrugation 
associated with frame bending of the cross 
sections (see eq. ( 11 )) 


strain energy per unit length of corrugation 
associated with longitudinal extension (see 
eq. ( 2 )) 


strain energy per unit length of corrugation 
associated with middle surface shear (see 
eq. (5)) 

strain energy per unit length of corrugation 
associated with torsion (see eq. (9)) 


longitudinal displacement 

one-half the relative shearing displacement of 
two adjacent trough lines (see fig. 3(b)) 

longitudinal displacement (function of z) along 
junction line (lj (see fig. 3(b)) 

longitudinal displacement (function of z) along 
junction line (2) (see fig. 3(b)) 


imaginary part of Rj, negative of imaginary part 
of R 2 when R-^ and R 2 are complex (see eqs. (B44) ) 

functions of z defined by equations (33) 


parameters defining the deformation of the cross 
section in its own plane (functions of z) (see 
fig. 3(c)) 

function of z defined by equation (C48) 

function of z defined by equation (D44) 

real part of R. and Rg when R 5 and Rg are complex 
(see eqs. (B44)) 



transverse coordinate (see fig. 1 (b)) 

imaginary part of R 5 , negative of imaginary part 
of R 5 when R 5 and Rg are complex (see eqs. (B44)) 

longitudinal coordinate (see fig. 1 (b)) 

defined by equations (B61) 

defined by equation (13a) 

defined by equations (B61) 

defined by equations (B61) 

defined by equations (B61) 

defined by equation (C 8 a) 

shear strain 

shear strain in plate elements 01, 12, 23 
respectively 

computed from equation (B23) and (B32) ; 
representable by equations (B46) through 
(B49) when R^, R 2 , .... Rg are complex 

B C D E 

complex conjugates of y^, y^, y^, y^ respectively 

B C D E 

complex conjugates of y^, y^, y,., y,. respectively 

defined by equations (C29) 
defined by equations (D24) 
longitudinal strain 

longitudinal strain in plate elements 01 , 12, 23 
respectively 
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defined by equations (B2) 


angle between sides of corrugation and horizontal 
(see fig. 1(a)) 

defined by equations (B14) through (B18) 

Poisson's ratio, taken as .3 for numerical work 
defined by equations (B40) 

defined by equation (C37) 

defined by equation (D32) 
cross-sectional normal stress 

cross-sectional normal stresses (functions of z) 
along junctions (l) and (5) respectively 

extreme-fiber bending stresses (functions of z) 
at junctions (T) and ( 2 ) respectively, resulting 
from frame bending of the cross sections 

middle-surface shear stress 

middle-surface shear stresses in plate elements 
01, 12, 23 respectively 

extreme-fiber shear stresses due to twisting of 
the plate elements 01, 12, 23 respectively 

rate of twist 

rate of twist of plate elements 01, 12, 23 
respectively 

factor in shear-stiffness equations (29) and 
(B63a); defined by equation (B63b) 

factor in equation (C42) for shear stiffness; 
defined by equation (C43) 

factor in equation (D38) for shear stiffness; 
defined by equation (D39) 

relative shearing stiffness, i.e. the ratio of 
shear stiffness of the actual corrugation to that 
of an identical corrugation with uniform middle- 
surface shear strain 



ANALYSIS 


The plate Is assumed to be composed of Infinitely many identical 
corrugations, all deforming in an identical way, and the analysis may 
therefore be based on a single corrugation. The form of the corrugation 
is shown in figure 1. A single corrugation is considered to be that portion 
between adjacent trough lines (lines labeled mn in fig. 1). The numbering 
system for the salient points of the cross section of a single corrugation 
is shown in figure 3(a). The notation for the geometry of the corrugation 
and coordinate systems is shown in figures 1 and 3(a). 

Considering a single corrugation, the shearing of this corrugation 
is imagined to be effected by a longitudinal shift of the trough line at 
station (o) through a distance uq in the positive-z direction and a 
similar shift of the trough line at station (5) in the negative-z direction, 
as shown in figure 3(b). Thus the total shearing displacement of one 
trough line with respect to the other is 2uq . The end-points of tne 
trough lines (points m and n in fig. 1) are moved only longitudinally. 
However, the rest of the points of a trough line are permitted to move 
both longitudinally and laterally, subject to certain constraints arising 
from the symmetry of the corrugation with respect to a vertical plane 
through the crest line (m'n' in fig. 1), the antisymmetry of the imposed 
displacements with respect to this plane, the requirement of continuity 
between adjacent corrugations, and the requirement that all corrugations 
deform identically. These considerations and requirements lead to the 
following constraints on the deformations of the trough lines: 

(a) The longitudinal strain is zero everywhere along a trough line. 

(b) The vertical displacement is zero everywhere along a trough line. 

(c) The trough lines at stations (o) and (5) curve into identical 
shapes. 

Certain deductions can also be made regarding the mutual internal 
reaction acting along the common trough line between two adjacent 
corrugations. These lead to the conclusion that, while there are an 
unknown longitudinal shear flow distribution and an unknown vertical shear 
distribution along the trough lines at stations (5) and , there is no 
horizontal running tension nor any bending moment transferred from one 
corrugation to the next across a trough line. 

Figure 4 shows schematically the type of linkage system that can be 
imagined to exist along the edges of an isolated corrugation in order that 
the isolated corrugation satisfy the above conditions and represent a 
single one of the infinitely many corrugations of the corrugated plate. 

Assumption regarding longitudinal displacements . - The longitudinal 
(z-wise) displacements at stations (0) and (3) of any cross section are 
+UQ and -Uq , as already discussed. The longitudinal displacements of 
the other middle-surface points of the cross section are assumed to 
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vary linearly between stations. These longitudinal displacements are 
shown in figure 3(b), which also shows their assumed antisymmetrical 
nature consistent with the antisymmetrical nature of the presecribed 
displacements at stations ( 5 ) and (5) . Therefore the longitudinal 
displacements of all middle-surface points of the corrugation are defined 
by one prescribed displacement parameter ug and two unknown functions 
of z : u-^(z) and U 2 (z) . If the resultant longitudinal shearing force 

F (see fig. 1(b)) is regard as prescribed instead of ug , the latter 
will become an additional unknown. 

Assumptions regarding displacements in the plane of the cross 
section . - Especially near the ends, the cross sections can be expected 
to undergo significant flexural deformations in their own planes, somewhat 
in the manner of a rigid-jointed frame. Therefore the deformations of a 
cross section in its own plane will be assumed to be inextensional, as is 
done in frame analysis . Certain degrees of freedom will be assumed for 
the displacements of stations (l) through (5) , consistent with the 
expected antisymmetry of the deformation pattern, and the displacements 
between stations will be assumed identical with those of the corresponding 
rigid- jointed frame, hinged at stations (o) and ( 3 ) . 

Three degrees of freedom are sufficient for this purpose, and 
figure 3(c) shows the three selected, as viewed from the positive end 
of the z-axis. The first two of these degrees of freedom are the same 
as employed in reference 1. The third is a rigid-body translation, of 
amount vg(z) , required because in the present analysis the trough lines 
are permitted to curve in the horizontal (xz) plane. Thus the displacements 
in the plane of the cross section are completely defined by three unknown 
functions of z : v-^(z), V£(z) and Vg(z) . 

Middle-surface extensional strains . - Referring to the foregoing 
assumptions regarding longitudinal displacements, and using the 
coordinate system of figure 3(a), the longitudinal displacements u 
for all points of the middle surface can be expressed in terms of ug , 
u^(z) and U 2 (z) . The expressions for these longitudinal displacements 
are given in the second column of table 1. The corresponding extensional 
strains e are obtained by differentiating these displacements with 
respect to z , and the resulting expressions are given in the last column 
of table 1. Because the longitudinal strains are antisymmetrical with 
respect to the crest line, it suffices to consider explicitly only the 
three plate elements listed in table 1. 
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Middle^surface shear strains . - The shear strains y of the middle 
surface of the plate elements of the corrugation, obtained from both the 
longitudinal displacements and the displacements in the plane of the 
cross section, are given in table 2. It is seen that they are constant 
across the width of any plate element (i.e., independent of sj, 82 > and 
S 3 ) , as a result of the assumptions that u varies linearly between 
stations and the cross-sectional deformation is inextensional. 


TABLE 2. - SHEAR STRAINS 


Plate element 

Shear strain, y 

01 

"1 ' "o dv o _ . 

e dz ~ Y 1 

12 

k - dz (v O COS0) + dz (v l sin0) = y 2 

23 

“ f dz + d2 (V in0COS0) + dz (v 2 Sin0) = Y 3 


Rate of twist of the plate elements . - If the corrugation length 
(2b) is several times the pitch (p) , it can be argued that the torsional 
strain energy will be a small fraction of the strain energy due to the 
flexural (frame-like) deformations of the cross sections.* Therefore, 
in computing the torsional strain energy of a plate element, it is 
probably sufficiently accurate to assume a constant rate of twist across 
the width of the element rather than to consider the detailed variation of 
rate of twist across the width. This constant rate of twist will be taken 
as the overall rate of twist corresponding to the displacements of the 
longitudinal edges of the plate element. For example, the rate of twist 

of the plate element 01 will be taken as I , in accordance with 

the edge displacements shown in figure 3(c) for this plate element. The 
rates of twist <j> obtained in this manner are given in tabxe 3. 


*The numerical results of reference 1 are consistent with this deduction. 
The same deduction is arrived at in reference 2. 
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TABLE 3 . - RATES OF TWIST 


Plate element 


Rate 

of twist, 

01 

til 

dz(e 

L ) 5 +1 


12 

d 

v 2) d 

'V-.COS0-, 

dzi 

w k J " dzi 

k H2 

23 

d 

^v^sin^G- 

d r 2v 2 cos0 ) 

dz 

l f 

+ dz( f 1 " *3 


Strain energy components . - As in reference 1 , the strain energy of 
the corrugation is assumed to arise from the following four sources: 

(a) middle-surface extension of the plate elements in the longitudinal 
direction, (b) middle-surface shear of the plate elements, (c) twisting 
of the plate elements, and (d) frame- like bending of the cross sections. 
Expressions are developed below for the density (i.e. strain energy per 
unit length of corrugation) due to each of these sources. In developing 
these expressions, use is made of the fact that plate elements 34 and 45 
contribute the same strain energy as plate elements 12 and 01 respectively. 

(a) Strain energy due to middle-surface extension: The strain 

energy, per unit length of corrugation, due to the longitudinal strains of 
the middle surface is 


U 


ext 



2 . 
£ i ds i + 


k 

'o 


2 

e 2 ds 2 



2 , 

c 3 ds 3 


( 1 ) 


where E' is the Young's- modulus associated with longitudinal extension, 
Cy, E2j £3 are defined in table 1, and S3, S2, S3 are coordinate shown 
in figure 3 (a). In writing equation ( 1 ), middle surface normal stresses 
in the transverse direction have been assumed negligible. The prime on 
the Young’s modulus symbol is a tracer to distinguish this Young's modulus 
from the Young's modulus associated with frame bending. 


Substituting the expressions for £3, £3, and £3 from table 1 and 
performing the integrations, one obtains 


du, \2 


U ext = b ll 


' aU l' 

jr , 


du du f 

+ 2b- „ + b o0 •, 

12 dz dz 22[dz 


du 2 .) 2 


( 2 ) 
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where 


>U ■ f E ' te(1+ ?> 


>i 2 - ?*’<* 


b 22 = -g- E't( 2 k + f) 


> ( 3 ) 


(b) Strain energy due to middle-surface shear: The strain energy, 

per unit length of corrugation, due to the middle-surface shear strain is 


U sh = Gt( Y J e + + 7 Y^f) 


( 4 ) 


where G is the shear modulus associated with middle-surface shear and 

Yl» Y2* Y3 are the plate-element shear strains given in table 2 . Substituting 

for the strains their expressions from table 2, one obtains 

“sh * °00 u 0 + c ll"l + C 22 U 2 + 2c Ol"o U l + 2c 12 u l"2 


dv dv 

+ d_„u_ -=-2- + d„ .u. 


dv. 


+ d„_u„ 


00 0 dz 10 1 dz 20 2 dz 


dv. dv. dv. 

+ d^u, - — + d 01 u. — — + d,.u_ — 
11 1 dz 21 2 dz 22 2 dz 




'00^ dz 


'22[ dz 


dv. dv dv„ dv dv 1 dv. 

ze 10 dz dz e 20 dz dz * 12 dz dz 


( 5 ) 


where 


'00 


'22 


'01 


* G 


G(f + 2|) 


- G — 
e 


'll 


'12 


c 7 (1 + f> 




> (6a) 
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- 2Gt(l - cos0) 
d 22 = - 2Gt sinG 

d^ * 2Gt sin0(l - cos0) 


00 

2Gt 

d io = 

20 

2Gt(l - cos0) 

d ll " 


\ 


> ( 6 b) 


e* Q = Gt(e + k cos 2 0 + |-) 
e*i * Gt sln 2 0(k + y cos 2 0) 

e 22 = J Gtf sin2e 

e*^ = - — Gt(2k + f)sin0cos0 
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(c) Strain energy due to twisting of the plate elements: The strain 
energy, per unit length of corrugation, due to twisting of the plate elements 
is 


U 


tw 


G ' (J l' f> l 


+ J. 



(7) 


where G' is the shear modulus associated with 'torsion (the prime being 
a tracer to distinguish it from the shear modulus G associated with 
middle-surface shear); J^, J 2 > >13 are the torsion constants of plate 
elements 01, 12 and 23 respectively, considered as bars of narrow rectangular 
cross section; and (fi^, <f> 2 > *)>3 are the rates of twist given in table 3. 
Substituting the expressions from table 3 for the rates of twist and 


J 


1 





for the torsion constants, one obtains 


U 

tw 


( dv, \ 2 

( dv \ 

1 . 7 

2 

l[dz J 22l 

[dz 


+ 2e 


dV ! dv 2 
12 dz dz 


( 8 ) 


(9) 


17 


in which 
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'12 
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Vk 


2J 
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(d) Strain energy due to frame bending of the cross sections: 
Considering a unit length of corrugation to be a rigid- jointed frame 
whose joint displacements are a superposition of the three modes shown 
in figure 3(c), while the joints are permitted to rotate freely, one 
obtains for the strain energy a quadratic expression in and • 

(vq is absent from this expression because it represents a rigid body 
translation. ) 

The derivation of this expression is given in appendix A of reference 
1 and will not be repeated here. In that derivation a parameter a is 
used which has the value 0 or 1 according to whether joints (5) and (J) 
are hinged or clamped. Only the case a = 0 is pertinent to the present 
analysis. Setting a equal to zero gives the following expression for 
the strain energy, per unit length of corrugation, due to frame-like 
bending of the cross sections: 


U b = a ll v l + 2a 
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' V 1 V 2 + a 22 V 2 
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where 
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with 


3 = 12 |(1 + |) + 6 f (3 + 4 f ) (13a) 

A 11 “ 216f(f + 2f)(2f-f 3f) + 432(f) 2 (f + 2f)2 

*22 ‘ 216 t(f +6 f>' 2 t +3 f> +144( t> 2 [ 3( t )2 + 30 tf +45( 7> 2 ] 

+ 864 <t> 3 7<l + 2 f> 
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432 tf (2 t +3 f )+144( t )2 f (12 t +21 f> 
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A 12 " - 432 £«f +3 ?>< 2 t +3 f> 

- 86 4 (f) 2 Cf + 3 f)(f + 2 f) 

A 13 " 432 tf (2 t +3 f> + 864( f> 2 f ( t +2 f ) 
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- 1296 f f<2 § + 3 f) - 288(J) 2 f (15 f + 27 j) 

- 1728 <t> 3 M + 2 f> 


> (13b) 




The symbol D appearing in equations (12) represents the frame 
flexural stiffness of the corrugation per unit width of frame, i.e. per 
unit length of corrugation. A detailed discussion of the symbol D has 
been given in reference 1, where it is concluded that an appropriate 
value for D is the plate flexural stiffness, i.e. 


D 


Et' 


12(1 - v*) 


(14) 


where E and v are Young's modulus and Poisson's ratio, respectively. 
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Total strain energy . - The total strain energy D of a single 
corrugation can be obtained by Integrating the sum of the foregoing 
strain energy densities over the entire length of the corrugation. I.e., 


U 


rb 

-b 


(U 


ext 


+ U . + U 
sh 


tw 


U fe )dz 


(15) 


where U ext , U 8 h» Utw» ^b are defined by equations (2), (5), (9) and (11) 
respectively . 

Total potential energy . - The potential energy of the prescribed 
shearing forces F along the sides of a corrugation, whose relative 
longitudinal displacement is 2 uq , equals - 2 Fuq . Adding this to the 
above strain energy U gives the total potential energy (TPE) of a single 
corrugation as 


TPE = - 2Fu q -HJ 


(16) 


Minimization of the TPE . - The TPE as defined by equation (16) is a 
functional of uq> u^(z), U2(z), vq(z), v^(z), V2(z) . In accordance 
with the method of minimum total potential energy (ref. 3) the "best" 
values of these quantities will be those which minimize the TPE. To 
obtain these "best" values, the technique of variational calculus may 
be used to form the first variation of TPE due to variations in uq, u^(z), 
U2(z), vq(z), vi (z), V2(z) and equate it to zero. This will lead to a 
system of field equations (primarily differential equations) and boundary 
conditions defining uq, u^(z), ..., V 2 (z) . 


The detailed execution of this procedure is given in appendix A. 
The resulting field equations, equations (A12) and (All) of appendix A, 
are repeated here for convenience: 
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and 


2c oo u o b + 


'01 


f 


u^dz = 


( 18 ) 


Different sets of boundary conditions are obtained in appendix A, 
depending on the nature of the end attachments. If there are point 
attachments at the ends of trough lines only (fig. 2(a)) the boundary 
conditions at z = ±b are 
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dz 
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dV 0 dv l dV 2 

2e — — + 2e — - + 2e — — + d u =0 

20 dz 12 dz 22 dz 22 2 




If the attachments at the ends of the trough lines are wide as idealized 
in figure 2(d), it is only necessary to replace the first of equations 
(21) by the condition 


v 



( 22 ) 


21 



Finally, if the attachments are as shown in figure 2(b), namely point 
attachments at the ends of both the crest lines and the trough lines, 
the boundary conditions are 


du 


1 


dz 


0 



( 23 ) 


v 


0 


0 


(24) 


v^cos0 + v 2 = 0 


(25) 


dv 0 dv l dv 2 

2(e i0 ‘ e 2O COS0) dz — + 2(e ll " e 12 COS0) di" + 2(6 12 " e 22 COS0) dz _ 


+ d^Uj. + ( d 2i " d 22 cos0)u 2 = 0 (26) 


The physical meaning of these boundary conditions is discussed in appendix A. 

Solution of equations . - Essentially, the basic solution of the problem 
consists of solving equations (17) for u^(z), ^(z), vq(z), vi(z), V 2 (z) 
in terms of uq , subject to appropriate set of boundary conditions. The** 
solution for u^(z) is then substituted in equation (18), which then gives 
uq in terms of F or F in terms of uq • 

The system of simultaneous differential equations (17) is linear with 
constant coefficients, and it can therefore be solved in a straight-forward 
manner. The full details, of the solution are in appendix B, and only the 
main features of the solution (those needed for computational purpose) will 
be cited. here. 

The numerical realization of the solution requires that equation (B20) 
be solved for its eight non-zero roots, R]_, R 2 , . . . , Rg . In equation (B20) 
c is any characteristic length (c was taken equal to the pitch p in the 
subsequent calculations) and the coefficients kg, k2, etc. are functions of 
the ratios of the elastic constants and of the parameters defining the basic 
shape of the cross section, i.e. 0, t/c, e/c, k/c and f/c . These 
coefficients are defined by equations (B13) . Because only even powers of 
R appear in equation (B20) , four of these roots are the negatives of the 
other four, as stated in equations (B21a) . 
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For most; cross-sectional geometries of interest through Rg are 
complex, and therefore the further description of the solution will pertain 
to that case. For a description of the solution procedure in the general 
case the reader is referred to appendix B. 

When R^ through Rg are complex, they come in conjugate pairs, and 
therefore four of them have the representation shown in equations (B44) 
with U, V, X and Y real numbers. Equations (B44) and (B21a) provide a 
complete representation of all eight roots. 


With U, V, X and Y known, the displacement ratios u-^(z)/uq , etc. 
can be computed from equations (B50) in which 5^ and t , 2 are defined by 

B C D E 

equations (B2); an d ?2 are defined by equations (B40); P ’ * ' , 

qB,C,D,E ^ gB,C,D,E ^,B,C,D,E are obtained by solving equations (B23) 
for j=l and 5 and noting equations (B31) and (B46) through (B49); and 


A 1 /u q , A^/Ug , A,./Uq , a 7 /u q , E 1Q / u o are obtained by solving equations 
(B52), (B53) or (B54) , depending on the type of end attachments. 


Relationship between F and Uq. - Equations (B50) give the displacements 

due to a prescribed value of uq , i.e. a prescribed value of half the 
relative shearing displacement of the two sides of a corrugation. For 
determining the displacements produced by a prescribed shearing force F , 
one needs the relationship between F and uq . This relationship is given 
by equations (B63) when Rq through Rg are complex. The symbols sv , sy, 
etc. in equations (B63b) are defined by equations (B60) . 

Equation (B63a) gives the overall shearing stiffness F/2uq of a 
single corrugation. One can define a dimensionless shearing stiffness 
parameter fi as the ratio of the actual shearing stiffness of a single 
corrugation to that of an identical corrugation having continuous end 
attachment producing uniform shear strain throughout the corrugation. The 
uniform shear strain of the latter corrugation due to the relative shearing 
displacement 2uq of its sides is 2uq divided by the developed width 
p' = 2e + 2k + f . The shear force F' required to maintain the relative 
shearing displacement 2uq is therefore Gt’2b*2u()/p' , which implies the 
following shear stiffness for the continuously attached corrugation: 


F’ 


2u 


0 


Gt*2b 

P’ 


(27) 


The relative shearing stiffness of the given corrugation is defined as 


ft 


F22 "o 

F72u 0 


(28) 
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Eliminating numerator and denominator of this expression via equations 
(B63a) and (27) respectively, one obtains 


« = f? - < i+ £ + H>* 


( 29 ) 


Stresses . - With the displacements u- L (z), U2(z), etc. determined, 
the stresses in the corrugation can also be obtained. Expressions for 
the various stresses will now be given, again restricted to the case in 
which the non-zero roots of the characteristic equation are complex. In 
order to avoid lengthy equations, the following short-hand notations will 
be employed: 'v 
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The longitudinal normal stresses a (J) and a/J) , along junctions ® 
® (see fig. 3(a)) are given by 


du^ du 2 

a ® = E ’ dT and a ( 2 ) = E ’ dT 


(34) 
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Substituting for u^ and U2 their expressions (eqs. (B50a) and (B50b)) gives 
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These stresses are positive if tensile. 


The middle-surface shear stresses in the plate elements making up the 
corrugation can be determined from the shear strains of table 2 and the 
displacement equations (B50) . These shear stresses will be denoted by 
T0l> Ti2 a nd T23 respectively. The following expressions are obtained for 
them: 
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t 23 C 

Gu_ 


" 2 fU 2 + (i/ + A 4 Q B )^ + (A 4 P B - A^®)^ 


U 5 S B + V B )^f + - V B >^f + «2 V 


+ V 2 sin0 + VjSinecose - V Q 


(36c) 


The positive sense of tq^ is shown in figure 5. The positive senses of 
the other two shearing stresses are the same as that of Tq^ . 

The extreme-fiber shear stresses due to torsion will be denoted by 
T 01’ t 12 an< * t 23 f° r P^ ate elements 01, 12 and 23 respectively. They 
can be obtained by multiplying the rates of twist in table 3 by G’t 
(ref. 4) and substituting equations (B50) for the displacements. The 
resulting expressions are 




— V 
e 1 


- £(v 2 + V^cose) 


(37a) 


(37b) 



2 -j(V^sin^0 - V 2 cos0) 


(37c) 


The positive sense of T 01 is shown in figure 5 as an example to indicate 
the positive senses of all three torsional shearing stresses . 

Due to frame-like deformation of the cross sections transverse bending 
moments are developed which vary linearly between stations and are zero at 
stations ( 5 ) and ( 5 ) . The transverse bending moment (per unit length of 
corrugation) at any station other than ( 5 ) and (5) is a linear function of 
vi(z) and V£(z) (vq playing no role because it represents a rigid-body 
translation of the cross section). The associated extreme-fiber bending stress, 
obtained by multiplying the bending moment by 6/t^ , will also be a linear 
function of v^ and V£ . The extreme-fiber transverse bending stresses at 
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stations (l) and ( 2 ) will /be denoted by and and will be 

considered positive if they are tensile in the lower fibers (see fig. 5). 
Expressions for these stresses in terms of and V2 can be obtained 
directly from equations (56) of reference 1 with the parameter a therein 
set equal to zero. The resulting expressions are 
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where v^/uq 
and the g_^ 


and v 2 /uq are now given by equations (B50c) and (B50d) , 
matrix elements are defined as follows: 
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Special cases . - The above results apply to the general case in which 
none of the dimensions e, f, and k is zero (fig. 6(a)). Analyses for 
the special cases f = 0 (fig. 6(b)) and e = 0 (fig. 6(c)) are contained 
in appendixes C and D respectively. In these appendixes only the end 
conditions of figure 2(a) needed consideration. Those of figure 2(b) are 
equivalent to continuous attachment when e or f approaches zero, 
because of the resulting complete suppression of the deformation of the 
end cross sections in their own planes. The end conditions of figure 2(d) 
are similarly equivalent to continuous attachment if f approaches zero, 
or to those of figure 2(a) if e approaches zero. 

The main results from appendixes C and D will now be cited. 

The special case f = 0 leads to a sixth degree (rather than a tenth 
degree) characteristic equation, (C19) , with six roots in the form (C23) 
and (C24). The numerical computation of the four non-zero roots can be 
done by means of the quadratic formula. With the roots determined, the 
relationship between the arbitrary constant (Aj , Dj , and E j ) can be 
obtained from equations (C28) and (C29). Equations (C39), with Ap, Xj, 
and Ag defined by (C40) , then give the displacements up(z), V 2 (z) and 
vg(z) . Equations (C42) and (C43) give the basic shearing stiffness, and 
equations (C45) through (C49) give the stresses. 

The special case e = 0 also leads to a sixth degree characteristic 
equation, (D14), which has six roots (R]_, ...» Rg) with the properties 
shown in equations (D18) and (D19). The four non-zero roots can be easily 
determined from the quadratic factor of equation (D14). With^the roots 
known, the relationship between the arbitrary constants (B- , D- and E-j) 
can be obtained from equations (D21) , (D23) and (D24) . Equations (D34) , 
with B^, B^j, and Bg defined by (D35) , then give the displacements U 2 (z), 
V 2 (z), and vq(z) . Equations (D38) and (D40) give the absolute shearing 
stiffness and the relative shearing stiffness respectively, and equations 
(D41) through (D45) give the stresses. 


NUMERICAL RESULTS AND DISCUSSION 

The foregoing analysis was used to determine numerical results on 
shear stiffness, stresses and deformations for selected cross-sectional 
geometries and end-attachment conditions. Poisson's ratio v was taken 
as 0.3, G was taken as E/[2(l + v)] , and no distinction was made 
between E and E' , or G and G' . In order to keep the number of 
computations within reasonable bounds, the numerical studies were limited 
to the case of the so-called symmetrical corrugation, that is the case in 
which the trough and the crest plate elements have equal width (2e = f ) . 

The numerical results were obtained by means of the equations discussed 
in the previous section. For determining the non-zero roots of the charac- 
teristic equation (B20), subroutine P0LRT of the IBM 360 Scientific 
Subroutine Package was employed. This subroutine can handle complex as 
well as real roots. The characteristic length c was taken equal to the 
pitch p . 
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For the sake of maximum generality the stiffness, stresses, and 
deformations are represented by dimensionless parameters'. For a given 
basic shape of cross-sectional centerline — i.e., for given values of 
h/p, e/p and f/p, there are two additional dimensionless parameters needed 
to define completely the geometry of the corrugation (except for its 
absolute size). The most obvious choices for these parameters are b/p 
(semi- length divided by pitch) and t/p (thickness divided by pitch). It 
was found, however, that the use of (b/p) • (t/p)3/2 is preferable to 
b/p alone as the choice for the first parameter; for then the dimensionless 
stiffness, stresses and deformations turn out to be relatively insensitive 
to the second parameter, t/p . A similar result was observed in the case 
of trough lines held straight (ref. 1), but in Jthat case bt/p2 was the 
significant parameter corresponding to (b/p) (t/p)3/2 0 f the present case. 

Shear stiffness . - Figures 7 through 9 give the basic numerical results 
for shear stiffness. The results are given in terms of the relative shear 
stiffness parameter ft , defined as the ratio of the absolute shear stiffness 
F/2uq of the actual corrugation to that of an identical corrugation with 
continuous end attachment producing a state of uniform shear (eq. (27)). 

To convert the relative shear stiffness ft to the absolute shear stiffness 
F/2 uq , it is only necessary to multiply !! by 2 Gtb/p’ , in accordance 
with equations (28) and (27) . That is 
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2u 
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Gt • 2b 
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(41) 


In 

for the 


3/2 

these figures ft is given as a function of (b/p)* (t/p) 
following range of cross-sectional geometries: 


h/p = .1, .2, .3, ,4, .5 

f/p(=2e/p) = .1, .2, .3, .5 

t/p = .005 and ,015 


Figure 7 is for the case of point attachments at the ends of the 
trough lines only (fig. 2(a)), figure 8 for the case of point attachments 
at the ends of both the crest lines and the trough lines (fig. 2(b)), and 
figure 9 for the case of wide attachments at the ends of the trough lines 
only (fig. 2(d)). 

Figures 10, 11 and 12 present the same kind of information as figure 
7, 8 and 9 respectively, but use log-log scales rather than semi-log 
scales in order to show more clearly the relationship between ft and 
(b/p) • (t/p)3/2 the regions of very low ft (close to zero) and very 
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high ft (close to unity). In figures 10, 11 and 12 the curves generally 
have a kink at ft = .5 . To the left of the kink each curve gives ft as 
a function of (b/p)*(t/p)3/2 . to r ight of the kink it gives 1 - ft 

as a function of (b/p>(t/p)3/2. 

The range of geometrical parameters covered in figures 10, 11 and 12 
is as follows: 


h/p = .2 and .4 

f/p(= 2e/p) = .1, .2, .3, .5 


for t/p = .005 and .015 ; and 


h/p = .2 and .4 

f/p(= 2e/p) = .2, .3, .5 


for t/p = .050 . 

The curves for t/p = .005 and .015 duplicate the information given 
in parts (b) and (d) of figures 7, 8 and 9. However the curves for 
t/p = .050 give information which is not contained in figures 7, 8 and 9. 

The closeness of the solid and non-solid curves in figures 7 through 
12 shows that ft is virtually a function of (b/p)*(t/p)3/2 alone, i.e. 
relatively insensitive to t/p , except in the region of very low 
(perhaps impractically low) values of (b/p)*(t/p)3/2 . 

Comparison of figures 7 and 8 (or 10 and 11) shows that a significant 
increase of shear stiffness results from having point attachments at the 
ends of the crest lines in addition to point attachments at the ends of 
the trough lines. (In the case of trough lines held straight (ref. 1) the 
increase of stiffness due to the additional set of attachments was much 
less significant.) 

Comparison of figures 7 and 9 (or 10 and 12) shows that a much larger 
increase of shear stiffness is obtained by changing from point attachments 
to wide attachments at the ends of the trough lines . This increase is 
also an upper limit to the increase that can be expected as a result 
of one-sided interference, like that shown in figure 3 of reference 5 
(also reproduced as figure 5 of ref. 1), between the troughs and the end 
member to which they are attached. 

As is to be expected, figures 7 to 12 show that, all other things 
remaining constant, an increase of h or f will lead to a reduction of 
the relative shear stiffness ft . Since increasing h or f implies 
increase of the developed width p' , equation (41) shows that the 
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absolute shear stiffness F/2ug will experience an even greater percentage 
reduction than the relative shear stiffness ft . 

It is of interest to compare the shear stiffness obtained in the 
present case with that obtained in the case of trough lines held straight 
(ref. 1). Figure 13, 14 and 15 present such a comparison for a particular 
shape of cross-sectional centerline (h/p = f/p = 2e/p = .2) and two 
values of thickness parameter: t/p = .005 and .020. In these figures ft 

is plotted as a function of the length-to-pitch ratio (2b/p). The dashed 
curves are for the case of trough lines held straight, the solid curves 
for trough lines permitted to curve. It is seen that the shear stiffness 
in the latter case is appreciably lower than in the former, except for the 
very short and very long corrugations (i.e. very small or very large values 
of 2b/p) and for the case of wide attachments at the ends of the trough 
lines (fig. 15). The lowering of the stiffness due to allowing the trough 
lines to curve is more pronounced for the thinner corrugation (t/p = .005) 
than for the thicker one (t/p = .020). 

The above-discussed differences in ft , between the case of trough' 
lines held straight and trough lines permitted to curve, suggest that 
analyses which make the simplifying assumption that the trough lines 
(along with all other generators) remain straight* (e.g., refs. 6, 7 
and 2), may be appreciably in error for some ranges of geometries if the 
plate does not actually have some external constraint which forces the 
trough lines to remain straight. 

Displacement and stress patterns for a particular geometry . - The 
manner in which the displacements and stresses vary along the length of 
the corrugation for one particular geometry is shown in figures 16, 17 
and 18, one figure for each of the three end-attachment conditions considered 
in the present analysis. The geometry is defined by the following numerical 
values : 


h/p = .4 

f/p = 2e/p = .4 

t/p = .015 

(b/p) (t/p)3/2 = .02 


These imply a length-to-pitch ratio, 2b/p , of approximately 21.8 and a 
6 value of 76°. Those quantities selected for plotting in figures 16, 


*The assumption of straightness of the generators is usually present 
implicitly as a by-product of the assumption of inextensional deformation 
of the middle surface. 
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17 and 18 are dimensionless measures of the longitudinal displacements 
u^ and U2; the displacements vq, and v£ in the plane of the cross 
section; the middle-surface shear stresses tqi» t 22 and T23» the extreme- 
fiber torsional shear stresses Tq^> xj_2 an< * 3^3 » longitudinal normal 
stresses CT (I) and a (2) at stations (T) and (2J ; and the trahsverse 
extreme-fiber stresses a ^ and a ^ at stations (l) and (z) due to frame- 
like bending of the cross sections. 

These figures show a number of things , the main one being that 
stresses and deformations are far from uniform along the length of the 
corrugation — i.e., "end effects" due to the discrete nature of the end 
attachments can penetrate an appreciable distance in toward the central 
portion of the corrugation. These figures also indicate that the major 
stresses are likely to be bending stresses associated with frame bending 
of the end cross sections. The peaking of the tqi shear stress near 
the end, to a value much higher than its average value, indicates that 
local buckling of the trough plate element due to shear near the end of 
the corrugation may need to be considered in the proportioning of the 
corrugation. The longitudinal normal stresses a (p) and 0 ( 2 ) » though 
smaller than the maximum normal stress due to flexure of the cross 
section, are much larger than the longitudinal stresses obtained in the 
case of trough lines held straight (ref. 1 ), and they therefore show that 
the assumption of inextensible generators would be less valid in the 
present case than it would be in the case of reference 1. Due to the 
longitudinal normal stresses a (T) » 0 (T) » 0 ( 3 ) = “ a ( 2 )» a © = 

- o (p) , and 0 ( 5 ) = 0 at junctions ( 0 ) through ( 5 ) , the crest and trough 

plate elements are in a state of bending in their own planes and the inclined 
plate elements may be in a state of combined compression and bending in 
their own planes. Figures 16, 17 and 18 show that these stresses reach 
peak values near (but not at) the ends. Local buckling due to them may 
also be a factor requiring consideration in the design of the corrugation. 

The graphs of vq, v]_ and V 2 show that the displacements in the plane of 
the cross section are an order of magnitude larger than the longitudinal 
displacements, and the nonlinearity of these graphs contradicts the frequently 
used simplifying assumption that the generators of the corrugation remain 
straight lines. 

Figure 17 provides a partial check on the correctness of the analysis 
and calculations. With f = 2e and point attachments at the ends of both 
the crest lines and the trough lines, an additional element of symmetry is 
introduced which, on physical grounds, should lead to the following 
characteristics for the stresses and displacements: 


T = T 

01 23 


T 1 = - t 1 T T = 0 

01 23 12 


© ■ °© 


u 2 ' u 0 - U 1 


v^cosG + ^2 = 


0 
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(The conditions = 0 and v^cose + V 2 = 0 express the fact that the 

inclined plate elements, of width k , undergo zero rotation at every 
cross section.) Examination of figure 17 shows that all of these charac- 
teristics except the first (tqi = T 23 ) are satisfied very well by the 
numerical results. The curves of tqi and T 23 , which should coincide, 
differ in ordinate by a few percent in the region z/b = .4 to 1.0 . 

This discrepancy may be due to round-off errors in the rathen lengthy 
calculation of T 23 . 

Maximum stresses . - Figures 16, 17 and 18 and similar results (not 
shown) for other geometries provided an indication of what kinds of 
stresses are significant and where their maximums occur. These maximums 
were then computed for a much larger range of geometries , and the results 
are presented in figures 19, 20 and 21, one figure for each of the three 
kinds of end conditions. These results were obtained for the case t/p = 
.015 ; however the dimensionless parameters used as ordinate and abscissa 
in each graph were so selected as to make the curves virtually independent 
of t/p . Figures 19, 20 and 21 may therefore be used for values of t/p 
other than .015. 

The range of geometries covered in figures 19 and 20 (point attachments 
at the ends of the trough lines only or at the ends of both the trough 
lines and the crest lines) is as follows: 


h/p = .1, .3, .5 

f/p = .1, .2, .3, .5 
3/2 

(b/p)(t/p) ' = .00035 to 2.5 


Figure 21 (wide attachments at the ends of the trough lines) covers the 
same range plus h/p = .2 . 

The stress maximums selected for plotting in figures 19 and 20 are 
the end values (at z = b) of the following stresses: (I) the extreme- 

fiber frame bending stress 0 (l) > (H) the middle-surface shear stress 

tqj , and (III) the magnitude of the maximum resultant extreme-fiber shear 
stress in plate element 01, i.e. | xqi | + l T 0ll • 

For the case of wide attachments at the ends of the trough lines, 
the frame bending stress a ^ was found to be generally larger than 

a (T) » ant * fig ure 21 therefore gives the end value of 0 (2) rat her than 

°(l) ‘ similarl y» l T 12 l + l T i 2 l can be lar 8 er than l T 0l l + 1*^1 ; 
figure 21 therefore gives the former sum instead of the latter. Because 
there was found to be very little twisting of plate element 01 in the 
case of wide attachments, the latter sum essentially equals |t | a l° ne * 
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An extensive survey of the peak values of the longitudinal normal . 
stresses CT (T) an< * a ( 2 ) wou ld have required much more calculation than 
the other maximum stresses, because the values of z/b at which the 
peak longitudinal stresses occur cannot be specified in advance. Therefore 
no graphical data are presented for them. The isolated numerical results 
in table 4 may, however, be of some interest. Table 4 gives, for selected 
geometries and end conditions , the peak values of the dimensionless 
longitudinal normal stress parameters and the location (z/b) of these 
peak values . 


COMPARISON WITH EXPERIMENT AND WITH ROTHWELL'S THEORY 

Rothwell in reference 2 cites some experimental data on shear stiffness 
of trapezoidally corrugated plates given by Horsfall in reference 8. The 
following test-specimen dimensions are given by Rothwell: h = .373" 

p = 3.55", 2e = f = .75", 6 = 20°, length = 2b = 18". The attachments 
consisted of 1/4-inch diameter bolts at the ends of both the trough lines 
and the crest lines. 

The results of the experiments, as given by Rothwell, are represented 
by the circles in figure 22. The lower curve in figure 22 is the theoretical 
prediction of the present theory, assuming point attachments at the ends of 
both the crest lines and the trough lines and assuming an isotropic material 
with Poisson's ratio of 0.3. The upper curve is the theoretical prediction 
given by Rothwell, based on his theory, which assumes inextensional 
deformation for the middle surface (thereby implying that the generators 
remain straight lines) but makes an approximate correction for the shear 
strain of the middle surface. 

As is to be expected, the present theory, with its more degrees of 
freedom, predicts lower shear stiffnesses than Rothwell' s, but it also 
predicts lower shear stiffnesses than those obtained experimentally. This 
may be due to the finite width of the bolt heads used in the end attachments 
or perhaps to interference between the deformation of the crest and trough 
plate elements and the member to which the attachment is made. There is 
not enough detail in reference 2 about the experiment to permit a more 
definite assessment of the cause of the discrepancy, and the original 
source, reference 8, is not available at the time of this writing. 


CONCLUDING REMARKS 

A theoretical analysis (based on the method of minimum potential 
energy) and numerical results have been presented for the elastic shearing 
of a trapezoidally corrugated plate with discontinuous attachments at the 
ends of the corrugations. The present work is an extension of previous 
work (ref. 1) in which the same problem was considered but with the trough 
lines assumed to be held straight. Thus the present work is more 
nearly applicable to a corrugated plate by itself, while the previous 
work was more pertinent to a corrugated plate fastened to a flat plate. 
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TABLE 4. - PEAK VALUES OF a/ 


End 

Condition 


h = f = 2e b f 

n n n ^ 


a CD p 
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Three different kinds of end-attachment conditions have been assumed 
in the present work, as in reference 1: (a) Point attachments at the ends 
of the trough lines only, (b) point attachments at the ends of the trough 
lines and the crest lines ; and (c) wide attachments at the ends of the 
trough lines only, the width of the attachment being the full width of the 
plate element at the trough. 

Based on the analysis, numerical results have been presented for the 
effective shearing stiffness and certain maximum stresses for a wide range 
of geometries. A knowledge of these quantities is felt to be of importance 
in the design and stress analysis of trapezoidally corrugated plates 
intended as shear webs. 

The numerical results confirm the by now well known fact that the 
absence of continuous attachment at the ends of the corrugations can 
cause a marked lowering of the effective shear stiffness, even for 
corrugation lengths many times larger than the pitch. This lowering of the 
shear stiffness results from the large frame-like flexural deformations of 
the cross sections permitted by the discontinuous end attachment. Of the 
three kinds of end conditions considered, point attachments at the ends of 
the trough lines gave the lowest shear stiffness, as was to be expected. 

A moderate increase in stiffness was obtained for the case of point 
attachments at the ends of the crest lines and the trough lines , but a 
very marked increase was obtained by having wide attachments at the ends 
of the trough lines only. The case of wide attachments represents an 
upper limit to the constraint provided by point attachments at the ends of 
the trough lines plus one-sided interference between the troughs and the 
member to which the troughs are attached. 

Because of the discreteness of the end attachments, the stresses can 
be quite non-uniform along the length of the corrugation. In particular, 
the middle-surface shear stress in the trough plate element was observed 
to peak significantly at the ends of the corrugation. TJje most significant 
stress, from the point of view of magnitude, was found to be an extreme- 
fiber bending stress associated with the flexural deformation of the end 
cross section. However, data have also been presented for the maximum 
middle-surface shear stress and the maximum extreme-fiber shear stress 
(combination of middle-surface and torsional shear stress). A limited 
amount of numerical data (table 4) was presented on the maximum longitudinal 
normal stresses. The longitudinal normal stresses vanish at the ends but 
in the interior reach peak values which may exceed the maximum shear stress. 

Por a given basic shape of cross-sectional centerline, two additional 
dimensionless parameters are required to completely define the geometry of 
the corrugations to within a scale factor, e.g. a thickness parameter t/p 
and a length parameter 2b /p . The numerical work revealed, however, that 
if a certain combination of length and thickness were used as one of the 
parameters, then the dimensionless shear stiffness and dimensionless 
stresses would be virtually independent of the second parameter. The 
combination parameter that serves this purpose was found to be (b/p) (t/p)3/2. 
(The analogous parameter when trough lines are held straight (ref. 1) was 

bt/p2 .) 
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Except - for the very small or very large length- tb-pitch ratios and 
except for the case of wide attachments a significant difference In shear 
stiffness was found to exist between the case of trough lines held straight 
(ref. 1) and the present case of trough lines which are permitted to cu*.ve. 

This suggests that the assumption frequently made (or Implied) in the 
shearing analysis of corrugated plates, that the generators remain straight 
lines, may be questional le in some cases and for some range of geometries. 

In reference 1, it was found that the analytical results taking into 
account the torsional stiffness of the plate elements making up the corrugation 
were only very slightly different from those obtained neglecting the torsional 
stiffness. Calculations, the results of which have not been presented, show 
tnat the same phenomenon is true in the present case, thus tending to justify 
the simplified manner in which the torsional strain energy was included in 
the total potential energy expression — i.e. by the use of an average rate 
of twist across the width of each plate element rather than the detailed 
pointwise rate of twist. 

Inasmuch as the present analysis is based on the method of minimum 
potential energy, one could claim that it over estimates the shear stiffness 
were it not for the approximate treatment of the torsional strain energy 
and the assumed absence of interaction between the frame bending moments 
and the longitudinal curvatures. Because of these simplifications one can 
only claim that the shear stiffness is probably over-estimated. 

It would appear that worth-while avenues of future analytical work on 
the shearing of corrugated plates should include the extension of the 
present approach to (a) the case of one-sided interference between the 
corrugation ends and the member to which the ends are attached and (b) the 
case of curvilinear (e.g., circular arc) corrugation. The shearing of 
the circular-arc corrugation has been studied by McKenzie (ref. 6), but on 
the basis of the assumption that the generators remain straight and 
inextensible. As already noted, such an assumption may not be appropriate 
in all cases. 
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APPENDIX A 


VARIATION OF THE TPE 


Equation (16) in expanded form is 


TPE = 
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The first variation of the TPE due to the small variations 6 uq, 

Su i ( z ), 6u£ (z) , 6vq(z), 6vi(z), Sv£(z) in the displacements u 0> u^Cz), 
.... V2<z) is 


6 (TPE) - 2 


rb 

J 

f b 

d(6u 1 ) 


d(6u 2 ) 

-r 1 

d(6u 1 ) ”j 

“b 

^11 dz 

dz 

+b 

b 12| dz 

dz 

dz | 
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+ b 22 dz dz~ ^ dz 
rb 


+ 2 
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Where the derivative of a variation appears in the integrand of 
equation (A3) , integration by parts will transform such a term so that 
the integrand involves the variation itself, rather than the derivative 
of the variation, and will also introduce boundary terms. Using the 
first term as an example. 


[b du x 

d(6u ) 

■*r J — — 

du 

b U A 2 

rb a 

J-b dz 

— az - 
dz 

dz 6U 1 

2 

, J -b dz 

— D 


6u.dz 


Reducing all integrands in this fashion wherever possible, and utilizing 
the boundary condition 


Vq (±b) = 0 , 6 Vq (±b) = 0 i 

which applies to all three types of end conditions shown in figure 2, 
one obtains 
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(equation continued on next page) 
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(A5) 


Equation (A 5) Is valid as it stands for the case In which there are 
point attachments at the ends of the trough lines only (fig. 2(a)). If 
the attachment is one of the other two types shown in figure 2 , equation 
(A3) must be modified to take into account the implied constraints on 
6 v^ and 6 v 2 at z - ±b . 

Thus, if there are wide attachments at the ends of the trough lines 
(fig. 2 (d)), it follows that 


(<Sv ) = 0 

z=±b 


(A 6 ) 


and the next-to-the last term in equation (A5) must therefore be omitted. 
If there are point attachments at the ends of the crest lines as well as 
the trough lines (fig. 2 (b)), the resulting constraint against horizontal 
displacement of the crest attachment point is expressible as 


[(v.sin0)cos0 + v„sin0 - v_] = 0 (A7) 

1 2 0 z=±b 

Taking into account equation (A4) and considering that sin0 4 0 , this 
becomes 


(VjCosO + v 2 ) z=±b = 0 

or, in variational form. 


(A 8 ) 


(fiv^cos© + 6 v 2 ) * 0 

A z=±b 


(A9) 
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Thus, for the type of end attachment shown in figure 2(b), 6 v£ in the 
last term of equation (A5) should be replaced by - 6 vicos 8 . The last 
two terms of equation (A5) can then be combined, leading to the following 
form of the equation for 6 (TPE) : 
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(A10) 


in which the symbol (...) has been used to represent terms that are 
identical to the corresponding terms in expression (A5). 

Differential equations . - In order for the TPE to be a minimum, 

<5 (TPE) must vanish for all possible values of Suq 6u2, 6vq, 6v^, 

6 v 2 consistent with the constraints. Thus the coefficients of the latter 
five quantities in the various integrands of equation (A5) must individually 
vanish, as well as the entire integral coefficient of 6uq . This leads to 
the following system of one integral equation and five differential 
equations, which apply regardless of the type of end attachment: 
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Boundary conditions . - The vanishing of 6 (TPE) also requires that 
the boundary terms of the 6 (TPE) expression vanish identically for all 
variations in v]_, V2» u^ and U£ at z * ±b consistent with the 
constraints. Referring to equation (A5) , which applies to the case of 
point attachments at the ends of the trough lines only, it is seen that 
this requirement leads to the following boundary conditions: 
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at z = ±b . Inasmuch as b 11^22 “ ^12 does not vanish, the first two 
of these equations may be replaced by 
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at z = ±b . Thus, the complete set of boundary conditions corresponding 
to the case of point attachments at the ends of the trough lines only 
(fig. 2(a)) consists of equations (A13) , (A14) and 


v Q (±b) * 0 (A15) 

which is the first of equations (A4) . 

For the case of wide attachments at the ends of the trough lines 
(fig. 2(d)), expression (A5) applies but with the next-to-the-last term 
excluded. As a result, the first of equations (A13) is non-existent, 
and the condition 


v x (±b) = 0 (A16) 


is use'd in its stead. Otherwise the boundary conditions are the same as 
for the previous case. 

For the third case, in which there are point attachments at the 
ends of the crest lines as well as the trough lines (fig. 2(b)), equation 
(A8) constitutes one of the boundary conditions, and equation (A15) a 
second. The remaining three boundary conditions, implied by the vanishing 
of the boundary terms of expression (A10) , are equations (A14) and the 
following: 


dv 0 dv i dv 2 

d ll U l + d 21 u 2 + 2e io dT + 2e ll d T + 2e i2 dT 


~ ( d 22 U 


dv C 

2 + 2e 20 dz 


+ 2e 


12 


dv i d M 

dz + 2 ®22 dz j 


icosO = 0 


(A17) 


at z = ±b . 

Equations (A14) in conjunction with the fact that dug/dz = 0 , are 
readily interpreted to mean that the longitudinal normal stress acting at 
the corrugation ends vanishes . 

Following the procedure of appendix B of reference 1, equations (A13) 
and (A17) can be shown to be equivalent to the requirement that certain 
effective in-plane shears at the ends of the plate elements vanish. 
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appendix b 


SOLUTION OF THE EQUATIONS FOR BASIC UNKNOWNS 


Equations (17) and (18) are the field equations of the physical 
problem. In this appendix equations (17) will be solved for ui(z), 
U 2 (z), vq(z), v^(z) and V 2 (z) in terms of uo , subject to boundary 
conditions of equations (19),. (20) and (21); or equations (19), (20). 
(22) and the second of equations (21); or equations (23), (24), (25) 
and (26) . Then the relationship between the shearing force F and 
the total relative shearing displacement 2uq of one side of the 
corrugation with respect to the other can be obtained from equation 
(18). Physical arguments can be used to show that u^ and U 2 should 
be even functions of z and vq, vj. and V 2 odd functions of z . The 
subsequent work will be simplified by considering only that solution 
of equations (17) which satisfies these conditions. 

Particular integral . - A particular integral of equations (17) 
will be sought in the following form, consistent with the even-ness 
of u^ and U 2 and the oddness of vq, v^, V 2 : 


u^ = constant , U£ = constant , Vq = 0 , v^ = 0 , V 2 = 0 

For this form of particular integral the last three of equations (17) 
are identically satisfied while the first two reduce to 


whence 


with 


C 11 U 1 + C 12 U 2 = " C Q1 U 0 


C 12 U 1 + C 22 U 2 = ° 


U 1 ? l u 0 


u 2 C 2 U 0 


— c c 
01 22 


C 11 C 22 " C 12 


C 01 C 12 


C 11 C 22 " C 12 


\ 

* (Bl) 

j 


> (B2) 


46 


J 



Thus a particular solution of equations (17) is 


u 


1 


? 1 U 0 



(B3) 


C haracteristic equation for complementary solutions . - The complete 
solution of equations (17) consists of two parts: (i) the particular 

integral, equations (B3) , and (ii) the complementary solutions, i.e. all 
the linearly independent solutions of the homogeneous system obtained by 
setting uq equal to zero . Solutions of this homogeneous system will 
first be sought in the following form: 


u 


u 


v 


v 


V 


1 

2 

1 

2 

0 


Ae 


rz 


Be 


rz 


Ce 


rz 


De 


rz 


Ee 


rz 


(B4) 


Substituting these assumptions in equations (17) with uq set equal to 


zero leads 

to the following restrictions on A, 

B * C , D , 

, 2 
b ll r - C ll 

, 2 
b r — c 
12 12 

2 d ll r 

0 

1 j 

2 d 10 

b 12 r2 - c 12 

, 2 

b 22 r C 22 

~ 2 d 2l r 

” 2 d 22 r 

1 , 

2 d 20 

- i d n r 

- 2 d 21 r 

a lf e ll r2 

a 12" e l2 r2 

2 

" e io r 

0 

~ 2 d 22 r 

a 12" e l2 r2 

a 22 -e 22 r 

2 

" e 20 r 

~ 2 d 10 r 

“ 2 d 20 r 

2 

" e io r 

2 

e 20 r 

2 

" e oo r 



A 


0 


B 


0 


C 

= 

0 


D 


0 


E 


0 


(B5) 
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Oi> 0i> o*> o. > cu> Hi > mi o > a*> 


Thus, for non-trival solutions of the form of equations (B4) , r must 
satisfy the following characteristic equation: 


b U r2 - c ll 

b 12 r2 ' c 12 

- \ d ll r 

0 

-i* 

, 2 
b 12 r -°12 

. 2 
b 22 r _C 22 

‘ 2 d 21 r 

“ 2 d 22 r 

-■h 

“ 2 d ll r 

_ 2 d 21 r 

a ir e n r2 

a 12 _e 12 r 

- e io 

0 

“ 2 d 22 r 

a 12' e 12 r2 

a 22 -e 22 r 

■ e 20 

- 1 d 10 r 

“ 2 d 20 r 

2 

- e io r 

‘ e 20 r2 

" e oo : 


= 0 


(B6) 


Expanding the determinant by cofactors based upon the fifth column, 
multiplying through by -32 for convenience, and introducing the short- 
hand notation 


■ a 12 


OQ > 

in 

■ b L 


h = 

2 

' C 12 


(_J. > 

III 

2 

' e 12 


III 

< ^ 

“ b ll C 22 

b 22 C ll 

III 

< 0 

' b 12 d 20 

III 

vD 
< T) 

4a l2 b 22 

* C 22 d 10 

III 

< 

4a 22 b 22 


\ 


12 12 11 22 22 11 


12 22 22 21 


11 22 12 21 


21 22 22 12 


21 12 22 11 


'22 12 


'22 22 22 


3 " 8c 12 e 20 “ 2d 10 d 22 


4 " C ll d 20 ' C 22 d 10 


5 " e 12®20 ~ ®10 e 22 


d 8 “ a 22 e 10 a 12 e 20 


d 9 ” 2d 10 d 21 ' 2d ll d 20 ” 8c 12 e 10 


> (B7) 


> (B8) 


a/ 
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(B9) 


n l “ b 12 d 10 “ b ll d 20 


n 3 " 2d 10 d ll ” 8c ll 6 10 


n 2 “ C ll d 20 - c 12 d 10 


n 4 ~ a 22 e 10 " a 12 e 20 


P 1 " e ll e 20 " e 10 e 12 


p 5 " 4a 12 b 12 + 4c 12 e 12 " d ll d 22 


p 2 ~ 4a ll b 22 + 4c 22 e ll ’ d 2l 


P 6 " b 12 d 10 “ b ll d 20 


p 3 B 4a ll b 12 + 4c 12 e ll “ d ll d 21 P 7 " C 12 d 20 " C 22 d 


10 


p 4 " 4a 12 b 22 + 4c 22 e 12 _ d 22 d 21 p 8 " a 12 e l0 ' a ll e 20 


q l ~ 2b 12 C 12 b ll C 22 ' b 22 C ll q 3 _ b ll d 21 b 12 d ll 


q 2 " d 22 e 20 ~ d 20 e 22 


q 4 - d 22 e i0 " d 20 e i2 


one converts equation (B6) to 


(k Q + k £ r 2 + k^r 4 + k 6 r 6 + k g r 8 )r 2 = 0 


where 


k 0 = 2e 00 X 5 + d 10 X 9 _ d 20 A 13 

k 2 = 2e 00 A 4 + d 10 A 8 “ d 20 A 12 ' 2e 20 A 17 + 2e lO A 21 
k 4 = 2e 00 A 3 + d 10 A 7 " d 20 A ll " 2e 20 A 16 + 2e 10 A 20 
k 6 = 2e OO A 2 + d 10 A 6 ” d 20 X 10 ' 2e 20 A 15 + 2e 10 A 19 


k 8 2e 00 A l 


2e 20 A 14 + 2e 10 A 18 


(BIO) 


(Bll) 


(B12) 


(B13) 
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with 


= 16b e 


16 (bg + ef) + ^ k ^ b n d 21 ~ 2b 12 d ll^ “ 4b ll d 22 m 


+ 4b 22 d ii e 22 


X 3 = 16 (ab + ce + fg) + Ztb 11 ( d 2 x h " d 22^ + 4c ll^ d 22 m ” d 21 k ^ 


2 /,2 


+ 8d i;L (c 12 k - b 12 h) + d 11 < d 22 “ 4a 22 b 22 4c 22 e 22^ 


X 4 16(af + eg) + 4h(2c 12 d 1;L c 11 d 2 i^ + 4 ^ c H d 22^ + a 22 C 22 d ll^ 


= 16ac 


8ed l + 8b 12 (e 10 k " e 20 m) + 8b 22 d ll d 5 


A 7 = 8gd 1 + «ed 2 + 8b 12 (e 1() h - e^j) + md 3 + kd Q + 2d n (e 1() d 7 - 
A g = 8^ + 8gd 4 + Jd 3 + hd g - 8d n c 22 d 8 


A„ = 8ad, 


10 


8e ni + 8b 11 (e 10 k - e^m) + Sb^d^ 


'll 


8en 2 + 8gn 1 + 8b n (e 1() h - e 2Q j) + kn 3 + Sc^^m 


+ 8d ll b 12 n 4 + 8d ll C 12 (e 22 e 10 " e 12 e 20 ) 


12 


8an i + 8 8 n 2 + 8c ll e 20 j + hn 3 ‘ 8c 12 d ll n 4 + 2a 22 d ll d 20 


A is 8an 2 




> (B14) 


)> (B15) 


> (B16) 
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X 14 " 16bp l 


X 15 4mp 6 + 16 ^ b 12 c 12 " °U b 22^ p l * 4b ll 6 20 P 2 

+ 4h 12 e 20 p 3 - 4d 11 e 12 d 1 + 4d 11 e 2 0 (b 22 d ll - b 12 d 21 ) 

+ 4< “lO (b ll p 4 ' b 12 p 5> 

A A A /b A A A 

X 16 * 4mn 2 + 4jp 6 + 4e 20 (c ll p 2 ’ C 12 P 3 } + 16(b 12 C 12 " b ll C 22 )p 8 


+ (d u d 10 - Ac n e 10 )p 4 + (4c 12 e 10 - d n d 20 )p 5 


+ 4d ll e 20 (c 12 d 21 c 22 d ll ) 


> (B17) 


A 17 = 4jn 2 + 16cp 8 + 4a 12 d ll P 7 


X 18 " 16bd 5 


A 19 = 16bd 8 + 16d 5 q l + 4q 2 q 3 “ 4b ll d 22 q 4 


+ 4b 12 d 10 (d 21 e 22 " d 22 e 12 ) ~ 4b 22 d 10 d ll e 22 


A 20 = 16d 5 C + 16d 8 q l + 4(G 12 d ll ' C ll d 21 )q 2 

+ 4a 22 d 20 q 3 + 4c ll d 22 q 4 + 4b 12 d 10 h 


+ d 10 d ll d 7 4a 12 b ll d 20 d 22 


> (B18) 


X 21 = 16cn 4 + 4hd 4 + 4a 22 d ll p 7 




51 



From equation (B4) it is evident that the quantity r appearing as 
the unknown in the characteristic equation (B12) has the dimension of 
(length) - ^. For computational purposes one should convert the characteristic 
equation to an alternate form in which the unknown is dimensionless. This can 
be done by introducing any characteristic length c (e.g., c could be 
taken equal to the pitch p ) and defining a dimensionless variable R as 
follows : 


cr 


(B19) 


In terms of R equation (B12) becomes 


2 ^44 

(k Q + R + -f R + 

c c c 


^ e 6 + !| r V - 0 


(B20) 


The roots of equation (B12) will be denoted by nj_, r2» . r^g, those 
of equation (B20) by Rj_, R£, . .., R^o, and these two sets of roots have 
th«> following relationship, in accordance with equation (B19): 


Rj = cr^ (j = 1,2,..., 10) (B19 ' ) 

The nature of the roots of the characteristic equation . - Because 
equation (B20) contains only even powers of R and has R? as a factor, 
the following properties can be postulated for the Rj : 


R 


4 



R„ 



R 


3 




(B21a) 


R 


9 


= 0 


R 


10 


= 0 


(B21b) 


Examination of equation (B20) shows that for a given cross-sectional 
shape (i.e. fixed values of 9, e/c, k/c, and f/c ) and fixed ratios of 
the elastic constants to each other, the Rj are functions of t/c only. 

The non-zero roots, Ri through Rg , may be real or complex, and for 
most geometries of interest they are generally complex. In that case the 
following additional properties may be ascribed to the Rj : 


R 


2 


V 


R 6 = V 


where R^* denotes the complex conjugate of R^ . 


(B22) 
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Complementary solutions associated with R^, R2, .... Rg. - The 

values of A, B, C, D, and G in equations (B4) associated with the 

particular non-zero root r ■ rj (j - 1, 2 8) will be denoted by 

Aj , Bj , Cj , Dj , and Ej . The relationships existing among these 
five coefficients can be obtained by substituting r = rj into equations 
(B5) and solving the last four of these equations for the ratios of B, C 
D and E to A . These equations can be written as follows: 


L 11 

L l ? 

L 13 

L 

L 

L „ 

12 

22 

23 

L 

L 

L 

13 

23 

33 

L 14 

L 24 

L 34 


where 
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* r 2 G * 
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E 1 

22 R j E 22 
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1 


L 12 


~ 2 

E d 21 R j 

r 


1 


L 13 


~ 2 

E d 22 R j 

r 


1 

G * 

L 14 


' 2 

E d 20 R j 



* 1 


L 22 


a ll' 

IcJ " [_E 6 11 

L 23 

= 

41 

t^2 T g ** 

,c] " [_E e l2 

L 24 

- 

G 

E 

** r, 

e 10 R j 

L 33 

B 

a 22 1 

'tl2 Tg ** 

wcj _ |_ E 622 

L 34 


G 

E 

** R 2 

e 20 R J 

L 44 

- 

G 
“ E 

e** R 2 
00 j 



+ §HrJ 2 *u>I 

+ f (ff 


£ * _ El * 2 
E 12 E b 12j 


1 £ d * R 

2 E d ll R j 


1 £ d * R 

2 E d 10 R j 


(B23) 
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with 


11 


12(1 - or 


<C/e) 2_2 k + A 22(t) cos2e + 4A 33(f) 2 sin40 


- ^ cos6 + 2 A 22 ^ Y sin 2 6cos0 - 2A^ 3 ^ sin 2 0 


’] 


12 


(c/e)~ 


12(1 - v 2 )g 2 


Mtf 


cos0 - ' sin 9cos0 


• I A 12 t ' A 23 t f( cos2e - sin2e > + A 13 f cos6 


22 


(c/e)' 


12(1 - v 2 )g 2 


[- 22 (tJ : 


+ 4 A 33 (f) cos 2 0 - 2A 23 f f cose 


■] 


‘u * M 


12 


22 


lk 
6 c 


S ♦ f) 


> (B26) 


J 


'11 


'12 


'22 


10 

* 

*20 


£ + £ 
e k 


c_ 

k 


£ + 2 - 
k f 


* * 
d ll “ d 22 


- 2 sin6 


d 21 = 2 sin0 (1 - cos0) 


-2(1 - cos0) 
2(1 - cose) 


> (327) 


> (B28) 


> (B25) 
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** e k 2. 1 £ 

= — + — COS 0 + -T ~ 

00 c c 2 c 


AA 

> 

'11 


, 2jk .If 2_ 

= sin Q[- + j - cos 0j 


** 1 f . 2. 

e 22 " Ic sin6 


A* 
5 10 


” f ( 2 ^ + ~) S±n0COS0 


AA 

*20 


-i -- sin6 
2 c 


aa If 2 

— sin 0cos0 
12 2 c 


'11 


- Ife 

3le 


c 2 c 4 

^ cos 6 + 2 ^ sin 0 


'22 


|(f + 2 f cos 2 ej 


'12 


- ■ 2 f sin2e ) cos9 




> (B29) 


> (B30) 
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The solution of equations (B23) will be denoted as follows: 


V A j 


1 

C j /A d 


c 

Y j 

V A J 

n 

D 

Y j 

v\ 


E 

hi 


(B31) 


From the last four of equations (B5) , it is seen that if r is replaced 
by its negative, B/A remains unchanged while C/A , D/A and E/A merely 
change sign. From this property and equations (B21a) it follows that 
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B 


Y 


4 


B 


Y 


1 


C 


Y 


4 



D 


Y 


4 



E 


Y 


4 



B 


Y 


3 


B 


Y 


2 





B 


Y 


8 


B 


Y 


5 






A 


> (B32) 




By summing the eight solutions of the form of equations (B4) 
corresponding to r - r^ through r = rg , expressing each exponential 
in terms of hyperbolic functions, discarding those terms not having 
the proper even-ness or oddness in z , and taking into account 
equations (B32), one arrives at the following part of the complementary 
solution: 

R i z * R_z * R S Z * R fi z 

u. = A* cosh + A„ cosh + A c cosh + A., cosh 

11 c3 c5 c / c 


u 


2 


* B R l z 

A 1 Y 1 COSh ~ 


, B R 2 Z 

+ A 3 y 2 cosh — — 


* B R 5 z * B R 6 z 

+ A c y,. cosh + A^y, cosh 

5 5 c 7 6 c 


* C R 1 Z * C R 2 Z * C R 5 z * C R 6 Z 

v, = A,y, sinh h A_Yo sinh + A c Yc sinh + A y, siiih 

1 11 c 32 c 5 5 c 76 c 


* D 


R l z , * D 


R 2 Z _ * D 


v 2 = Vl Slnh ~ + A 3 Y 2 Sinh ~ + A 5 Y 5 Slnh ~ + A 7 Y 6 


R 5 Z . .* D R 6 Z 

sinh 


* E 


R 1 Z . .* E 


R 2 Z , .* E 


R^z 


* E . R 6 Z 


v. = A-y, sinh + A„Yt sinh + A,y c sinh + A n y, sinh 

0 l'l c 3 2 c 5 ' 5 c 7 6 c j 


where A^, A^, A^, A^ are certain linear combinations of the 
may be regarded as new arbitrary constants. 


A '< 

j * 


and 


Complementary solution associated with the root R = Rg = 0. - 

Substitution of r = rg = 0 into equations (B5) gives the following 
relations which the coefficients Ag, Bg, ...» Eg must satisfy: 


> (B33) 
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__ «. 


(B 34 ) 


2 9 

Since c n c 22 “ c 12 an< ^ a ll a 22 “ a \2 
have the solution 


are not zero, these equations 


b 9 = c 9 - d 9 = 0 


Eg = indeterminate 


\ 


(B 35 ) 


Substitution of these values, together with r = Vg = 0 , into equations 
(B 4 ) gives the following complementary solution: 


u„ 


= v„ 


= 0 


v 


0 


E 

9 


\ 

► (B 36 ) 


where E„ 

9 


is an arbitrary constant. 


This solution gives Vq even in z rather than odd. The constant 
Eg may therefore be equated to zero. Thus the root R = R9 = 0 makes no 
constribution to the complementary solutions. 


Ad ditional comple mentary solution not of the form of equations (B 4 ) . - 
The presence of a repeated root (see eqs . (B 21 b)) of the characteristic 
equation indicates that there exists a complementary solution that is not 
of the form of equations (B 4 ). This complementary solution can be found by 
inspection if one assumes the following form for it: 


u, = A 


10 


u„ = B 


10 


c io z/c 

D io z/c 


E io z/c 


(B 37 ) 
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where Aig, B 10» •••* E 10 are constants. Equations (B37) satisfy the 
symmetry and antisymmetry requirements with respect to z . Substituting 
equations (B37) into the original differential equations (17) (with the 
uq term omitted) leads to the following conditions on A^g, B 10» •••» B 10 
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(B38) 


whence 


where 




c „ = 

D „ = 0 

10 

10 

E 10 = 

indeterminate 

A 

E„E „ 

10 

nio 

B 

E E 

10 

^2 10 


> (B39) 




C 12 d 20 ~ C 22 d 10 
2C(C 11 C 22 ~ C 12 ) 


C 12 d 10 C ll d 2Q 
2C(C 11 C 22 ~ C 12 ) 


(B40) 


58 



Thus the complementary solution of the form of equations (B37) is 


u 2 “ 5 2 E 10 

- 0 


E io 2/c 


(B41) 


where is an arbitrary constant. 

Complete solution . - The complete solution of the differential 
equations (17) (satisfying the symmetry and antisymmetry requirements), 
obtained by adding together the particular integral (B3) and the 
complementary solutions (B33) and (B41) , is 


c. u. + A, cosh + A„ cosh + A r cosh 

101 c3 c5 c 


* R fi z 

+ A ? cosh — + 


* B R 1 Z * B R 2 Z * B R 5 Z 

U 2 * S U 0 + A 1 Y 1 COSh ~ c~ + A 3 Y 2 COsh T" + A 5 Y 5 cosl1 ~ c~ 


* B R 6 Z 

+ A ? y 6 cosh — + C 2 E 1q 


* C R 1 Z * C R 2 Z * C R 5 Z * C R 6 Z 

v 1 = A 1 y 1 sinh + A^ sinh — + A^ sinh — + A ? y 6 sinh 

* D R l z * D R 2 Z * D R 5 Z * D R 6 Z 

v 2 = A 1 y 1 sinh — - + A^ sinh — - + A^ sinh + A ? y 6 sinh 


* E R 1 Z * E R 2 Z * E R 5 Z 

A-V- sinh + A„y~ sinh + A r Y«- sinh 

11 c 3 2 c 5 5 c 


+ a*y E Sinh ~ + E 10 f 


> (B42) 
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* * * * 

The unknown constants A^, A 3 , Ag, A 7 and E^q are determined in 
terms of uo through the boundary conditions, which are: equations (19), 

( 20 ) and ( 21 ) for the case of point attachments at the ends of the trough 
lines only; (23) through (26) in the case of point attachments at the ends 
of the trough lines and the crest lines; and (19), (20), (22) and the 
second of equations ( 21 ) in the case of wide attachments at the ends of 
the trough l<nes. 


Substitution of the expression for uq(z) from equations (B42) into 
equation (18) gives 


<c 00 + C 01 5 l ,b + c oil 


K c (l ¥/3= . „ E 2 b + A 5 c ( . V 

— - ■=— sinh + — — sinh + — — sinh 

v u 0 R 1 c u 0 R 2 C u 0 R 5 C 


A_ R,b E „ > 

+ — £- sinh — + — C,b 
u 0 R 6 c u 0 1 


2 u„ 


(B43) 


With A-^/uq , etc. determined through the boundary conditions, equation 
(B43) yields the effective overall shearing stiffness F/2 uq of a single 
corrugation. 


Special form of solution when through Rg are complex. - The 

procedure described above is quite general; it applies regardless of 
whether the eight pon-zero roots of the characteristic equation (B20) 
are real or complex. For almost all cross-sectional geometries of 
interest, however, the eight non-zero roots of this equation turn out 
to be complex. It is therefore worth-while to investigate the special 
form taken on by equations (B42) and the boundary condition equations 
in that case. 


Considering the 
equations (B22), one 


case of complex roots, 
can represent R]_, R 2 , 


and taking into account 
R 5 and R 5 in the form 


R^ = U + iV R 5 = X + 1Y 

R. = U - iV R, = X - iY 

2 . D 


(B441 


where U, V, X and Y are real numbers. Furthermore, from the last four 
equations of (B5) it is evident that if r is replaced by its complex 
conjugate then B/A, C/A, D/A and E/A are also changed to their complex 
conjugates. Applying this to the complex conjugate pairs rj_, ^ and 
r 5 , rg , and taking into account equations (B31) , it follows that 
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B C D E 

Thus the y * * * which appear in equations (B42) have the following 
representation ’ 


y-, = 


P B + iQ B 


P B - iQ B 


S B + iT® 


S B - iT B 


> (B46) 


P C + iQ C 


S C + iT C 


P C - iQ C 


S C - iT C 


Yo - 


p D - iq D 


„D D 
S - iT 


P E + iQ E 


S E + iT E 


P E - iQ E 


S E - iT E 


B B 

where P , Q , etc, are real numbers. 


(B47) 


P D + iQ° 


S D + iT D 


> (B48) 


> (B49) 


Substituting expressions (B44) and (B46) through (B49) into 
equations (B42) gives the following form of the complete solution, 
applicable to the case in which the eight non-zero roots of the 
characteristic equation are complex: 


u l(z) A i Uz Vz . A 4 . , Uz . Vz 

— c, + — cosh — cos — + — sinh — sin — 

u 0 1 u 0 c c u 0 c c 


, A 5 . Xz Yz , A 8 . . Xz . Yz . E 10 _ 

+ — cosh — cos — + — sinh — sin — + 

u o c c u o c c u 0 1 


(B50a) 
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u 2 (z) 


? 2 + 


( A 1 , A 4 a b] . Uz Vz 

— - P + — Q | cosh — cos — 

U 0 U 0 J 


K B 

A > 

1 

, L Uz 

j Vz , 

r * 5 

B A 8 -b! 

+ — p 

Q 

sinh — 

sin — + 


S + — T 

l u o 

u o J 

c 

c 

.“o 

U 0 J 


( — S® - — T®1 

\ u o u o ) 


. , Xz . Yz ± 10 _ 

slnh — sir 1- 

c c u 0 2 


(B50b) 


v 1 (z) 


/ A 

4 

A i cl 

cosh 

Uz . 

Vz . 

( a 

A i _C 

A 

4 

— P 

=■ Q 

— sin 

— + 

— P 

+ — Q 

i u o 

u 0 J 


c 

c 

l u 0 

u o 


. , Uz Vz 
smh — cos — 
c c 


( A „ 

A > 


f A 

A 

8 C 

5 m C 

^ Xz . Yz , 

5 „C 

, 8 m C 

— S 

— T 

cosh — sm — + 

— S 

H — T 

K 

u o J 

c c 

l u o 

u o 


. , Xz Yz 
sinh — cos — 
c c 


(B50c) 


v 2 (z) 

'S P D 

- — Q° 

cosh 

Uz 

a ,_ Vz 

_L 

(1 

-ip° 

A N 

+ ^q d 

u o J 

u o 

l u o 

“o Q J 

c 

sm — 
c 

T 

i u o P 


. , Uz Vz 
sinh — cos — 
c c 


+ 

fk 

-8 s D 

A > 

5 _D 
— T 

,Xz . Yz , 
cosh — sm 1- 

(a 

s D + 1* T D ) 


. u o 

u o 

c c 

, u 0 

"0 1 


Xz Yz 

— cos — 
c c 


(B50d) 


v o (z) 


f — P E - — Q^| cos h — sin — + 

"d J 


^ U 0 


f> A ps + 1* 

u 0 u 0 


. , Uz Vz 
sinh — cos — 
c c 


(a 

8 E 

A > 

A 5 _E 

Xz . Yz , 

5 0 E 

A ^ 

8 _E 

— S 

l u 0 

T 

u o 

cosh — sin H 

c c | 

— S 

l u o 

+ — T 

u o 


. , Xz Yz 

sinh cos — 
c c 


10 z 

u o c 


(B50e) 
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where A^, A^, A 5 and Ag are new arbitrary constants related as follows 
to the previous ones: 


A 


1 


* 




* * 
A 5 + A 7 



A 8 



(B51) 


These four arbitrary constants and the fifth one, E^q , are to be 
evaluated from the boundary condition equations listed in the previous 
section after equations (B42) . These lead to the following sets of 


simultaneous equations for cosh 


A 4 , Ub A 5 , Xb 

— cosh — , — cosh — , 
u 0 C u 0 


, Xb 
cosh — 
c 

, E io 

and 

u o 

• 



N u 

N 12 

N 13 

N 14 

0 

N 21 

N 22 

N 23 

N 24 

0 

N 31 

N 32 

N 33 

N 34 

N 35 

N 41 

N 42 

N 43 

N 44 

N 45 

N 51 

N 52 

N 53 

N 54 

N 55 


— — 


— — 

A 1 , Ub 

— cosh — 


0 

=0 c 

A 4 , Ub 

— cosh — 


0 

u 0 

*5 . Xb 

— cosh — 


0 

"0 

A 8 u Xb 

— cosh — 


N / 

u o 

E 

10 

u o 


4 

N 5 

— — 


— ^ 


(B52) 


for the case of point attachments at the ends of the trough lines only 
(fig. 2(a)); 
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The matrix elements in these equations are defined as follows: 





N 


N l] 

ae 

U tu cv - V sv N 21 = pBN ll " qBn 12 



N 12 

- 

U sv + V tu cv N 22 = pBN 12 + ^11 



N 13 

- 

a /> B B 

X tx 9y - Y sy N 23 = S B N 13 - T B N 14 


>• (B55a) 

N 14 

= 

rv o B B 

X sy + Y tx cy = S*N 14 + T 

) 


N 31 

= 

_E „E r\ 

P tu cv - Q sv 



N 32 

= 

_E r\ , r\ ts 

P sv + Q tu cv 



N 33 


_E r\ r\ J o 

S tx cy - T sy 


> (B55b) 

N 34 

= 

E r\ m E o a 

S sy + T tx cy 



N 35 

= 

b/ c 

y 





> 


N 41 

= 

d ll^ + d 21 a l + e ll 6 l + e 12 S l + e 10 6 l 



N 42 

= 

~ r\ r\ ~B~C~D~E 
d ll tu sv + d 2l“2 + e ll g 2 + e 12 B 2 + e 10 B 2 



N 43 

= 

~ r\ ~B~C~D~E 
d ll Cy + d 2l“3 + e ll B 3 + e 12 S 3 + e 10 6 3 


> (B55c) 

N 44 

= 

d ll tX sy + d 21 a 4 + e ll g 4 + e 12 B 4 + e 10 B 4 



N 45 

= 

d ll C l + d 21 ? 2 + 6 10 

> 





> 


N 51 

= 

d 22 a l + e l2 B l + e 22 6 l + e 20 6 l 



N 52 

= 

d 22 a 2 + e 12 B 2 + e 22 B 2 + e 20 6 2 



N 53 


d 22 a 3 + e 12 P 3 + e 22 B 3 + e 20 B 3 


> (B55d) 


(equation continued on next page) 
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(B56) 


II 

m 

53 

d 22°4 + 

e B C + e 0° + e B E 
e 12 P 4 e 22 P 4 e 20 P 4 


N 55 = 

d 22 ? 2 + 

®20 

> 


N. 


C l d ll “ C 2 d 21 


N. 


? 2 d 22 




— ,_.C o ^ n C A. n . ,_D f\ ^ n D 

= (P tu cv - Q sv)cos6 + (P tu cv - Q sv) 


N ^2 = (P^ sv + tu cv)cos0 + (P^ sv + tu cv) 


= (S C tu cv - T C sc)cos0 + (S° tu cv - 1° sv) 


= (S C sv + T C tu cv)cos0 + (S° sv + Q° tu cv) 


} (B57a) 


N 51 

= 

*41 “ 

N^cosO 

¥ 52 

= 

N 42- 

N^cosO 

W 53 

= 

N 43 

N^^cos0 

*54 

= 

*44 " 

N c , cos0 
54 

*55 

s 

N 45- 

N,_,.cos0 

*5 

= 

N 4 ' N ! 

.COS0 

) 

*41 

s 

P tu 

r* _C 

cv - Q 

f 42 

* 

„C r\ 
P sv 

C ^ 
+ Q tu 

*43 

m 

e C r\ 
S tx 

r\ m C 
cy - T 


ae 

c C 

S sy 

m C rs 
+ T tx 


> (B57b) 


(B58) 


> (B59) 
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where 


*\ 


r\ 


J Vb 

sv 


sin — 
c 

n 


Yb 

sy 


sin — 
c 

A 


Vb 

CV 


COS 

c 

r\ 


Yb 

cy 

— 

COS 

c 



, Ub 

tu 

= 

tanh — 
c 

r\ 


, Xb 

tx 

E 

tanh — 
c 


> (B60) 


J 


and 


1 

B 

a 2 = 

B 

a 3 

B 

« 4 = 
= 

•l ■ 
- 

»; ■ 

•j - 

So = 


B r\ „B a n 
P cv - Q tu sv 


„B a n J n 
P tu sv + Q cv 


S B cy - T B tx sy 


r o , „,B r\ 
S tx sy + T cy 


U (P C cv - Q C tu sv) - V(P C tu sv + Q° cv) 


U (P C tu sv + Q C cv) + V (P C cv - Q C tu sv) 


X(S C cy - T C tx sy) - Y(S C tx sy + T C cy) 


X(S C tx sy + T C cy) + Y(S C cy - X C tx sy) 


U(P^ cv - Q C tu sv) - V(P D tu sv + Q D cv) 


U (P D tu sv + Q D cv) + V (P D cv - Q D tu sv) 




> (B61) 
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(equation continued on next page) 
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e° = X(S D cy - T D tx sy) - Y(S D tx sy + T° cy) 


S? = X(S° tx sy + T C cy) + Y(S° cy - T° tx sy) 


8® = U(P E ot - Q E tu sv) - V (P E tu sv + Q E cv) 


3^ = U(P E Ci £ + Q E Or) + V (P E ^ - Q E Cu&) 

8 E = X(S E cy - T E tx Q) - Y(S E tx sy + T E cy) 

g E = X(S E tx sy + T E cy) + Y(S E cy - T E tx sy) 


J 


and 


d Xl = d 22 = ' slne 


21 


'11 


'22 


'12 


'10 


sin0(l - cos©) 

+ tt + 
^ + 2 f cos2e ) 

if fife 


c 2 c 4 ^ 

f cos^e + 2 ^ sinV 
k f 


j f sin 2 e + 


> (B62) 


If 2 

~r — sin 0cos0 + 

2 c 


c 2 1 

cos© - 2 — sin 0cos0 

I J 


- i(2 — + — )sin0cos0 
2 c c 


'20 


1 f . A 
- y ^ sin0 


Substituting expressions (B44) into (B43) , and making use of (B51) and 
the definitions of Cqq and cqi (eqs. (6a)), one obtains the following 
formula for computing the shearing stiffness F/2uq of a single 
corrugation in the case of the non-zero roots being complex: 



(B63a) 


68 



where 


! _ r -30 g -£ 
C 1 u 0 «1 b 


(k l Ub 

— cosh — 

u 0 


„ r\ 


U tu cv 

__ 

+ V sv 

. 2 

„2 

U + 

V 


^4 Ub' 

U sv - V tu cv . 

(I > 

5 . Xb 

X tx cv + Y sy 

— cosh — 

u n c 

V 0 J 

2 2 + 

U z + V 

— cosn — 

L u o c 

2 2 
X + Y 


(A 

u, 


\ n n ^ 

8 , Xb X sy - Y tx cy 

— cosh r *- 

0 c X 2 + Y 2 


(B63b) 


It will be noted that in equations (B51) , (B52) and (B53) the 
^1 Ub 

combinations — cosh — , etc. are regarded as the unknowns rather than 
u 0 

A 1 

— , etc. alone. This is done in order to avoid having extremely large 


matrix elements in the simultaneous equations when b/c is large. 
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APPENDIX C 


SPECIAL CASE: f = 0 


In this appendix the special case f = 0 will be considered, in 
which junctions lines @ and (5) coincide and form a line of points 
of inflection of the cross sections (see fig. 6(b)). Only the end- 
attachment conditions of figure 2(a) (i.e. point attachments at the 
trough line ends only) need be considered in conjunction with the case 
f = 0 , for any of the other two types of end attachment would completely 
prevent cross-sectional deformation and therefore be tantamount to 
continuous attachment. 

Along the common junction line formed by junction ( 2 ) and (5) when 
f = 0 the longitudinal deformation must vanish, and the vertical 
displacement must also vanish. These constraints can be expressed as 
follows: 


u 2 (z) = 0 (Cl) 

and 

[v^ (z)sin0]sin0 - [v 2 (z)]cos0 = 0 (C2) 

The variational form of these equations is 

6u 2 = 0 (C3) 

o 

Sv^ = 6v 2 cos0/sin 0 (C4) 


Introducing the foregoing conditions into equation (A5) gives the 
following expression for the first variation of the TPE: 


5 (TPE) - (6u„) | (2c 00 u 0 + 2 c q1 u 1 - ^)dz 


-b 


+ 2 


Mb ^ 

-bC 11 H, 2 


d u , 1 dv Q 

, . + c 01 u 0 + C 11 U 1 + 2 d io dT 

dz 


i dv 

. 1 , 2 cos 

+ — d 


2 11 dz . 2 

sxn 


(Su^)dz 


- 2 


f b fl . ^1 

_ b [^2 10 dz 


+ e 


^ V 0 , , cos0 , V 2 

+ ( e m T~ + e on ) 


00 . 2 
dz 


'10 . 2 20 M 2 

sm 0 dz 


1 


(6v rt )dz 


(equation continued on next page) 
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cos0 


+ 2 I I (a,, C0 S t~ + 2a 


f \o 

' -b L 


+ a 00 )v. 


11 , 4 q 12 , Z D 22 2 

sin 0 sin 0 


1 j du l COS0 , COS0 . N d v 0 

- 2 d n IT XT 7 ' (e io XT' + e 20 ) 7T~ 


sin 0 


sin 0 


dz 


/ 2 \ 
, , cos_0 + 2e ' » cos ^ + ~ t « 1 

11 4 12 2 

k sin 0 sin 0 22/ 


?V 

dz 2 


(dv^idz 


du„ 


+ (2 b 


li ST Su i> 


-b 


COS0 , 

llUl sin 2 e I " 10 


{* 


cos0 


sin 2 e e20 J dz 


+ 2 (e ' ' + 2e” + e”)^ («v 2 ) 


11 sin 4 0 


'12 , 2 . ' ~22'dz 

sin 0 


-b 


(C5) 


where 


'00 


'10 


'20 


'11 


'22 


'12 


Gt(e + kcos 2 0 + ~) 


- — Gt(2k + f)sin0 cos0 


- — Gtf sin0 

o f 2 1 d 2 2 

Gt sin 0 (k + j cos 0) + G 1 — ^ j cos 0 

V e Z k 

1 2 d 2 
j Gtf sin 0 + G'-f 

k 

12 d 2 

-T- Gtf sin 0 cos0 + G' — cos0 

Z k 




> (C6) 


At this stage the vanishing of f has not yet been introduced 
into equation (C5). In order to incorporate this condition, f may be 
allowed to approach zero in all terms of equation (C5) except those 
associated with frame bending, that is an, a ^2 and a22 • The strain 
energy of frame bending for the case f = 0 cannot be obtained 
correctly by letting f -*• 0 in the equations for a^, a i2 an d a 22 — 
equations (12) and (13). The reason for this is that if condition 
(C2) is first imposed (as it has been) to prevent vertical displacements 
of stations (2) and ( 5 ) , the subsequent imposition of the condition 
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f -*■ 0 will lead to clamping (zero rotation) at the vertex formed by 
stations ( 5 ) and ( 5 ) as they meet, rather than to the condition of free 
rotation corresponding to the point of inflection (zero moment) which 
should be present at this vertex. In order to obtain correctly the 
zero moment condition existing at the vertex, f must be allowed to 
approach infinity, rather than zero, in those terms of equation (C5) 
which arise from strain energy of frame bending, namely an, ai2» 
and a 22 • Doing this, one obtains the following limiting values of an, 
a^i anc* a22 for use in equation (C5) : 



a 22 -.2 3 A 22 (kj " a 22 

Be 

where 


B 

5 12 1 (1 + 1> 



(C8a) 

A 11 

i 432(ff(l + 

I) 



A 22 

2 «2(f) 3 a + 



> (C8b) 

A l2 

5 - 864(f) 3 <l 


) 



Incorporating the above limiting values of a^, a^J^ and a 22 
into equation (C5) and letting f approach zero in equations (C6) 
leads to the following expression for 6 (TPE) : 


5 (TPE) 


(<S V 


-b 


(2c 00 U 0 + 2c 01 U l " b )dz 


+ 2 


d u. 


i . ts. 


(~ b ll 2 + C 01 U 0 + C 11 U 1 f 2 d 10 dz 

-n v dz 




(equation continued on next page) 
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- 2 




,2 a* 

d v 0 , cose d v 2 ,,, 

+ ®00 2 e 10 22 Vrt)dz 
dz sin 6 dz 


l J< 5 v 


f fo 


2 

. _ , , , cos .8 , cos 9 , 

+ 2 (a,, — + 2a 10 — + a oo> v 


11 . 4 * 12 , 2 . 22 ' 2 

sin 6 sin 8 


JL ,4 ^ U 1 cos 8 ~ , cos 8 d v 0 

" 2 11 dz , 2 ~ e 10 . 2 q . 2 
sin 8 sin 8 dz 


cos 8 . ", cos8 


d 2 v. 


- (e^ + 2e' ^==- +-ej 2 )“r K«v,)d« 

11 sin 8 X sin 8 dz 


} 


du l 

+ « b ll dz « U 1> 


-b 


[ d ll U l 


cos8 ~ , cos8 dv Q 

. 2. 10 . 2. dz 

sin 8 sin 8 


+ 2 1 e ' — — — + 2e' co - se + e' 1 — — 

Z|e ll , 4„ X 12 . 2 22jdz 

- sin 8 -I 


! ( ; i: 


2 q 
cos 8 


sin 8 


1 ( 6 v 2 ) 


-b 


where 


e^Q - Gt(e + k cos 8) 




- - Gtk sin8 cos8 


= Gtk sin 8 + G 


■& + S coe2e ) 


'22 


'12 


-S 

J 2 

G’ cose 
k 2 
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(C9) 


> (CIO) 


L 



Differential equations and boundary conditions . - The vanishing of 
6(TPE), equation (C9), provides the following equations governing u-^(z), 
v 0 (z), V£(z) and F : 


4 c oo u o b + 2c oi 


ll u i dz ' 


2F = 0 


(Gil) 


d 2 u 


- b 


1 . , , 1 , dv 0.1, cos0 dv 2 

+ C„,u„ + C^U, + 7T d n „ - — + — d, — 


11 . 2 '■010 11 1 2 10 dz 2 11 . 2 dz 


\ 

0 


dz 


2 2 
du d v ~ d v A 

2 “10 dz 00 , 2 ^10 , 2„ . 2 


sin 0 


id 


> (C12) 


dz 


sin 0 dz 


d d 2 

cos0 U 1 cos6 v 0 


d 2 V. 


a 22 V 2 2 d ll . 2 Q dz " 6 io . 2 q , 2 e 22 , 2 

sin 6 sm 0 dz dz 


= 0 


du^ 

dz 


= 0 


z=±b 


r dv dv 

, cos0 , „ , cos0 0 , „ , 2 

d, , s- u, + 2e' =— T— + 2e ' 


11 -r 2 fl '1 
sin 9 


'10 . 2 dz 

sin 0 


'22 dz 


= 0 


z=±b 


(C13a) 


where 


22 


2 

cos 9 , „ cos0 

a u ~rr * 2a i2 ~X 


+ a 


sin 9 


sin 0 


22 


e 22 6 


• S°±J + 2e’ cos6 + e’ 

11 , 4 Ze 12 . 2. 22 

sin 0 sm 0 


> (C14) 
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I 


Besides equations (C13a) , the boundary condition 


v Q (±b) = 0 


(Cl 3b) 


of the general case applies also to the present case. 


Particular solution of the differential equations . - A particular 
solution of equations (C12) having the proper parity with respect to z 
is 

C 01 

u. u 

1 c u o 


> (C15) 


Complementary solutions of exponential form . - The complementary 
solution, which is the solution of the homogeneous system obtained by 
omitting the term cq^uq from equations (C12) , is first assumed in the 
form 


u^ = A'e 


rz 


v 2 = D'e 


rz 


v Q = E’e 


rz 


> (C16) 


Substituting these expressions into equations (C12)_with the uq term 
excluded, one obtains the following conditions on A', D', E' and r : 


b 



c 


11 



11 


cos 8 
2 

sin 6 


I A 

~ 2 d ll 


cos 8 
sin^0 



“ 2 d 10 r 



cos 8 ~2 

, 2 . r 
sin 8 



d io r 



cos 8 
sin^0 




(C17) 
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1 


which lead to the following characteristic equation for r : 


b ll r ~ C 11 


1 , COS 0 

- 2 d ll T ^7 r 

sin 6 


1 , cos 6 ~ , , '2 

2 d ll 2 r a 22 " ®22 

sin 0 


I d 10 r 


cos0 ~2 

— e r 

1° . 2 e 
sin 0 


" 2 d 10 r 


cos0 ~2 

‘ e 10 , 2 r 
sin 0 


- ;• r 2 
00 


Equations (C18) when expanded becomes 



where 


20 


22 


40 


k 20 + k ?2 - 


( c ) 2 ]* 2 + ^40 + *42 (if 



= 0 


R = cr 


k 02 2 a 22|_ C ll e 00 


r * 


- Ko > 2 ] 


3 2 

G cos 8 . 

■ O 

E sxn 0 


*o * "1 ? 

( a io * ” c n e ooJ s;Ln e 


r * ^ ^ i 1 ** a * '*i 

l^ll^io^ 4 (d ll ) e 00 2 d ll d 10 e loJ 

" 7 e 2 2 [ C ll e 00 - i (d 10> ] 

EU'G) 2 * r~* k . 2 0 '* ,2~l cos 2 e 

EIE2 b llj e 00 c dln (e 10 ) . 4 

«- Jsin 0 


E’ G * 

E E b ll a 22 e 00 


E'fGl * 

k 42 EiEj b ll e 00 e 22 


= 0 (C18) 


(C19) 


(C20) 


> (C21) 
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with 




a 


22 


(c/e)' 


12(1 - 


2 '2 . 4. 
v )3 sin 0 


E 


A n cos 


e - 


A 12 k COS0 + A 22 


(tf] 




e oo 

n 

O | (T> 

, k 2 

■1 cos 0 

c 

e io 

= - 

— sin0 cos0 
c 

~* 

1 

G' 1 [ c 

e 22 

3 

E sin 4 0 ^ 6 


> (C22) 


Denoting the roots of equation (C19) by R^, R 2 » ...» Rg» the following 
properties may be postulated for them: 

R 2 = - R x R 4 = - R 3 (C23) 


R 5 = R fi = 0 (C24) 


The corresponding values of r will be denoted by r^, r 2 , . . . , rg . 

Thus, five solutions of the homogeneous system having the form of 
equations (C16) exist, corresponding to the five different roots of 
equation (C19). For each such solution the relationships which must 
exist among the coefficients A', D', and E* can be obtained by 
substituting the_particular value of r into equations (C17). Denoting 
by A! , D! , and E! the values of the coefficients associated with the 
particular root r J = r^ , the following relationships are implied by the 
last two of equations (C17) for j = 1,2,3 or 4 : 


.2 , , ~ 2 
a 22 C “ e 22 (cr j ) 

cos0 , ~ . 2 
' sin 2 e < “j ) 


d:/a: 

3 3 


1 , COS0 , ~ x 

2 d ll c . 2 <“j> 

sin 0 J 

e' cose (c ; ) 2 

10 . 2/ c V 

sin 0 J 

- *00 Cor j )2 


E ! / A ! 
J 3 


— d c ( cr ) 
2 10 


— 


— — 


— 


(C25) 
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Substituting for a£ 2 » e i0' e 22» d 10 and d ll their definitions from 

equations (C14) , (CIO), and (6b), one converts equation (C25) to 


where 


L 11 L 12 


L 12 L 22 


D!/A! 
J J 


E'/A! 
J J 


1 G cose j* ~ 
2E si„ 2 e 11 ^ 


7 I ^ 


(C26) 


11 


- (*] 


* _ e * R 2 ) - — — 

22 22 ” - 


G k cos 2 6 Z2 

E c . 2 0 R j 
sin 6 J 


G COS0 T*. ^2 

'12 - E . 2 ~ e *0 E j 

sin 9 J 


> (C27) 


J 22 


G . J2 

e* R 

E 00 j 


The solution of equations (C26) is 


d:/a: 

J 3 


~D 

Y : 

e:/a: 
i j 
_ _ 

ii 

~E 

Y i 


(C28) 


where 


^ ■ 7 ! (l 22 - £2 12> 

J J sin 0 


^ *= 7^l^l d !o- £ l2^ d ?l )/(£ ll £ 22-^ 

J J sin 9 


> (C29) 
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From equations (C29) and (C23) it is evident that 


~D ~D ~D ~D ~E ~E ~E ~E 
Y 2 “ - *1 . Y 4 » -Y 3 , Y 2 - -Y X . Y 4 ’ - Y 3 


(C30) 


For j * 5 (r - r^ - 0) equations (C17) give 


r i 


— 


“ — • 

i 

o 

o 

o 


K 


0 

11 


5 



0 *22 0 




0 

0 0 0 


*5 


0 

— _ 




— — 


(C31) 


whence 


- 0 , E 5 = indeterminate 


(C32) 


_ Summing the^four solutions of the form (C16) corresponding to 
r «= r^ through r 4 and making use of equations (C28) and (C30) , one obtains 


— — 


— 



— 

U 1 


1 

1 

1 

1 

V 2 


~D 

Y 1 

_D 

" Y 1 

_D 

Y 3 

_D 
- Y 3 

V 0 


"E 

Y 1 

~E 

- Y 1 

'E 

Y 3 

-<i 

w w 

1 



1 C 



z_ 

1 C 



z^ 
3 c 



(C33) 


Expressing the exponential functions in terms of hyperbolic functions and 
discarding the terms which do not have the proper symmetry (in the case of 
u 3 ) or antisymmetry (in the case of v 2 and vq) with respect to z , one 
converts equations (C33) to 
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U, 


— R l z — R 3 z 

A, cosh + A_ cosh 

1 c 3 c 




R z 



R~ z 

A 1 

sinh 

X 

c 


A 3 

sinh 

c 



R l z 

~E 


V 

X 1 

sinh 

c 

+ y 3 

A 3 

sinh 

c 


(C34) 


where A. and A^ are new arbitrary constants arising from, certain linear 
combinations of the previous ones. 


From equations (C32) and (C16) the following solution corresponding 


to r 


0 is obtained: 


E 5 


Since Vq m * constant is even in z , rather than odd, the constant 
Ec may be equated to zero. Thus the pertinent solution contributed by 
the root t = r^ = 0 is identically zero. 

Complementary solution not of exponential form . - The existence of 
a repeated root (see eq. (C24)) indicates that the homogeneous system has 
an additional solution that is not of the form of equations (C16) . This 
additional solution can be obtained by inspection if one postulates that 
it has the following form consistent with the evenness of u^ and the 
oddness of V£ and vq 


v 2 - Dg z/c (C35) 

v 0 - E' z/c 
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Substitution of this assumption into equations (C12) with the uq term 
omitted shows that equations (C35) are indeed a solution of the 
homogeneous system provided that 


°6 


0 


where 



5 = - 2cc 11 /d 10 


(C36) 


(C37) 


Thus the following additional solution of the homogeneous system is 
ob tained : 

U 1 " A 6 

v 2 » 0 (C38) 

v 0 " « A 6 Z/C 


Complete solution of the differential equations (C12) . - Summing 
the particular Integral (C15) and the various solutions (C34) and (C38) 
of the homogeneous system, one obtains the following complete solution 
of equations (C12) having the pertinent symmetry and antisymmetry 
iroperties : 


C 01 — ^1 Z — ^3 Z 

u„ + A, cosh + A„ cosh + A’ 

c i;L 0 1 c 3 c 6 


~D 


V 


~D v 


2 “ Yi A j_ sinh + Y 3 a 3 sinh 


R 3 Z 


(C39) 


■E - 


R 1 Z . 'E - 




0 * *1 A x sinh — + y 3 A 3 


A„ sinh 


3- + TV ■=• 


6 c 


Evaluation of the arbitrary constants . - The three arbitrary 

constants A^, A^ and A£ can be determined from the boundary conditions, 

equations (C13a) and (C13b). The following equations are obtained which 
can be solved simultaneously for the three constants: 
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Ai Rjb 

PH P 12 0 - sinh — 0 


A 3 V 

P 21 P 22 / P 23 Uq 8inh c ° 


P 31 P 32 P 33 A 6 /u 0 


P = R 
*11 1 




*§•!.(!)'> 


, ~E "j. cos0 

+ e *o TX 

sin 0 




f;Vk cos^i + E - 

L Y ^ C sin 2 0 G e 22UJ/ 

-~ E ~fr COS0 "|l . R 3^ 

+ V!o C ° th “ 


p * d * Spa. 6 + 2e* ? ^ s6 - 

r 33 a H 2 10 ^ 2 

J J " L sin 0 sin 0 


_cos0/sin 0 
11 1 + e/k 


(C41b) 
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Relationship between F and u^. - With u^(z) determined in 

terms of uq (the first of eqs. (C39)), equation (Cll) can be used to 
obtain the following relationship between F and u^ : 



(C42) 


(C43) 


A relative shear stiffness ft can be defined the same way as in the 
general case. Equation (29) of the main body of the paper applies with 
f set equal to zero and i/> replaced by $ 


Stresses . - The longitudinal normal stresses are identically zero 
at stations (0) , (2) , (3) and (5) . The non-zero longitudinal normal 
stresses a /*f\ along junction iQL) are obtained from the equation 


a Q 



(C44) 


Eliminating u^ through the first of equations (C39), one obtains the 
following expression for the dimensionless stress parameter O/^-sc/E'Uq : 


o 





sinh 



sinh 


R 3 z 


(C45) 


The dimensionless middle-surface shear 
table 2 and equations (C39)» are 


stresses , 


as obtained from 
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T 01 c 

Gu„ 


X 12 C 

Gu. 


e 1 + k/e u 


(f-‘M 


, R 1 Z 
cosh 

) c 


-f-- 

u oi e 


~p - Y R~z A' / 

Y, R~ Icosh — t — + — ( 5 

3 3 ) c u 0 ? 


> (C46) 


= _ c/k A lf~P ~ 

1 + e/k + u q ( Y 1 R 1 


c ~E ~ \ R l z 

cote - i y, R, cos0 cosh 

k '1 1 c 


^Y3 R 3 c °te - t - Y3 R 3 cose 


cosh 


V 


-^(t + E cosa J 


J 


where the subscripts 01 and 12 stand for the plate elements 01 and 12 
respectively. 


From the rates of twist in table 3 and the displacement expressions 
(C39) the following expressions are obtained for the extreme-fiber shearing 
stresses due to the twisting of plate elements 01 and 12 respectively: 


T 01 C 

G'u 


t COS0 


e . 2 
sin 0 


W(z) 


T 12 C 

G'u 


t_ 

k 


• 2 
sm 0 


W(z) 


where 


A 1 ~D ~ ^l z A 3 ~D ~ R 3 Z 

W(z) = — y R cosh H — y_ R_cosh 

Uq 1 1 c Uq 3 3 c 


(C47) 


(C48) 


The frame bending moments and associated extreme fiber bending 
stresses are zero at stations (6) , (T) , (3) and (T) . Referring to 

equations (D44) of reference 1, one obtains the following dimensionless 
expression for the extreme fiber transverse bending stresses 
along junction (l) : ^ 
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(C49) 


Eu o 


'c 


1 - V 


t_ c_ 
k e 


8sin^0 


[■ 


£ cos 6 + 3 



in which a ^ ' is positive for compression in the upper fibers , 
Jn the lower fibers (see fig. 5). 
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APPENDIX D 


SPECIAL CASE: e = 0 

Figure 6(c) shows the special case e = 0 , in which the plate 
elements at the troughs are of zero width, with the result that the 
two adjacent inclined plate elements meet to form a vertex along the 
trough line. 

This special case can be obtained from the general case (appendix 
A) by first imposing along the junction line (T) the same displacement 
conditions as exist along the junction line (0) , namely 


u l(z) = u 0 

CD1) 

v 1 (z) = 0 

(D2) 


and introducing these conditions (and their variational forms 6u-^ = 

6uq, 6v^ = 0 ) into the 6(TPE) expression, equation (A5). This 
will eliminate from the 6(TPE) those terms arising from longitudinal 
extension and twisting of plate elements 01 and 45. 

In order to eliminate the terms associated with middle-surface 
shearing of these plate elements, the condition e -*■ 0 should then be 
introduced into all the remaining coefficients except those associated 
with the strain energy of frame bending (a^, a ]_2’ a 22^' 

The condition e -*■ 0 will not lead to the correct strain energy of 
frame bending because this condition, imposed after condition (D2), 
implies clamping (zero rotation) rather than free rotation along the 
trough lines in the limiting case. In order to obtain correctly the 
zero moment condition existing along the trough lines the condition 
e -*■ 00 must be imposed instead in those terms (a-Q, a 12> a 22^ associated 
with frame bending of the cross sections. 

Applying the above procedure to equation (A5) gives the following 
expression for 6(TPE): 


6 (TPE) = («u 0 ) (2c 00 u 0 + 2c 12 u 2 - F )dz 

— b 


+ 2 


— b ' 


d 2 u 2 ! dv 0 

+ c_ «U_ + C„„U„ + TT d„ 


, . 22 , 2 12 0 22 2 2 20 dz 

-b ' dz 


i d V 
+ 2 d 22 dT, 


(6u 2 )dz 

(equation continued on next page) 
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f b 

[ 1 

du 2 

c 

> 

d V 2 \ 

u 

(2 d 20 

dz 

+ a — 

10 . 2 

dz 

620 dz 2 J 


d v 2 \ 

a? - J (6v 0 )dz 


+ 2 


r b f- 1 . 

f^2_ 

d v 0 

d V 2 Y 

J__ b [ a 22 V 2 2 d 22 

dz 

' 20 dz 2 

622 dz 2 J 


[ 2b 22 ^< 6u 2>]|_ b 

[( d 22 u 2 + 2e 20 3T + 2 *22 ar)^]!^ 


(D3) 


where 


c 00 “ Gt/k 


(D4) 


‘00 


Gt (k cos 2 0 + j) 


(D5) 


22 


2(1 - v ) 


(fi 


3 |(1 + 2 | cos 6) 2 

mi 


(D6) 


Differential equations and boundary conditions . - From the vanishing 
of S(TPE), equation (D3) , the following equations governing U 2 , vq and 
V£ are obtained: 


4c 00 U 0 b + 2c 12 


-b 


u 2 dz - 2F = 0 


(D7) 


d 2 n 2 x dv 0 1 dv 2 

" b 22 , 2 + C 12 U 0 + C 22 U 2 + 2 d 20 dz + 2 d 22 dz 
dz 


= 0 


1 . du 2 
tt a 


d 2 v„ 


a 2 v„ 


+ 6, 


+ e„ 


2 20 dz "00 , 2 20 . 2 

dz dz 


> (D8) 


du„ 


d 2 v„ 


d 2 v„ 


„V„ - — d. 


*22 v 2 2 “22 dz " e 20 . 2 "22 . 2 

dz dz 


- e. 
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( du„ 

£ 

Jz=±b 


= 0 


( d 22 U 


dv 0 dv 2 
2 + 2e 20 dz + 2e 22 dz 


^=±b 


(D9a) 


= 0 


Besides equations (D9a) , the boundary condition 


Vq (±b) = 0 


(D9b) 


of the general case applies also to the present special case. 

The above development is for the case of point attachments at the ends 
of rhe trough lines only (fig. 2(a) or 2(d)). When e = 0 the presence of 
additional attachments at the ends of the crest lines (fig. 2(b)) would 
completely prevent cross-sectional deformation and therefore be tantamount 
to continuous attachment. 

Particular solution of the differential equations . - A particular 
solution of equation (D8) having the proper parity with respect to z is 


C 12 

c 22 “° 


v„ = 0 


= 0 


> (DIO) 


Complementary solutions of exponential form . - The complementary 
solution, which is the solution of the homogeneous system obtained by 
omitting the term c^ug f ron i equations (D8) , is first assumed in 
the form 


u 


v 


v 


2 

2 

0 


Be 


rz 


De 


rz 


Ee 


rz 




> (Dll) 


J 
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Substitution of these expressions into equations (D8) with the c^^O 
term omitted leads to the following conditions on B, D, £ and t : 


b 22 r 1 2 - c 22 

i , , 

" 2 d 22 r 

” 2 d 20 r 


A 

B 


0 

" 2 d 22 r 

a 22 ” e 22 r 

-2 

- e 20 r 


D 

- 

0 

1 

~ 2 d 20 r 

*2 

~ e 20 r 

-2 

" e oo r 


E 


0 


from which arises the following characteristic equation for r : 


b 22 r * C 22 


1 j " 

2 d 22 r 


1 j " 

2 d 20 r 


2 d 22 r 


- ~2 
a 22 “ e 22 r 


-e 20 r 


Equation (D13) when expanded becomes 


- 2 d 20 c 


e 20 r 


-2 

e oo r 


[Wlf] + [ £ 20 + £ 22(i) 2 J £2 + [ £ 40 + £ " 


where 


R = cr 


(D12) 


(D13) 


(D14) 


(D15) 


89 


02 



k 


22 


k 


40 


(I ) 2 - J <d !o )2 ] 

(IJ 3 [l (d 20 )2e 22 - J d 20 d ll e 20 + d !2 (e ?5 )2 ' ^‘38*!! 

+ i (d *i )2 *ssj 


(!) 2 ‘*22 - C $ 2 *$ 2 S SS ] 

rdf b h^mi - (e *o )2] 


— - I* b* e** 
E E 22 22 00 


(D16) 


42 


IL £ £l u* - - 

E E E b 22 00 22 


with 


a* 

22 


d ) 2 f (1 + 2 t dos9)2 

2(1 - v 2 ) (1 + 2 |) 


-** k 2 1 £ 

e** = — COS 6 + -r- — 

00 c 2 c 


(D17) 


Denoting the roots of equation (D14) by R-^, R 2 > •••, R 5 , the following 
properties may be postulated for them: 


R 2 " - R 1 


R 4 = - R 3 


(D18) 


R c = R, = 0 

5 6 


(D19) 


The corresponding values of r will be denoted by r-p , ...» r^ . 
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Thus, there exist five solutions of the homogeneous system having 
the form of equations (Dll) . For each such solution the relationships 
which must exist among the coefficients^ B, D, and E can be obtained 
by substituting the specific value of r into equations (D12) . Letting 
Bo, Da, and Eo denote the values of 6, D, and E associated with the 

^ j A ^ 0.1 £.11 I 1 -1.J 1__f J 1 J 3 1 A.L. — *1 A. 4-.«a 


- 2 
a 22 C 


- e 


22 V 

^2 


20 


<C V 


4 : 






A A 

Vj 


1 d 22 c(c V 




m. 


CM 

<H^ 

o 


A A 

y B J 


I d 20 c(cf j ) 


(D20) 


Substituting for &22> e 22» e 20» ^00» ^22 anc * ^20 t^eir definitions from 
equations (D5) , (D6) , (6b) and (A2), one converts equations (D20) to 


£ 11 

L 12 


V B j 


— — 
IP 

2 E d *l R j 

L 12 

^22 


y\ 

— 

1 r * 

2 E d 20 R j 


vD21) 


where 


"11 




"12 


- — e**R 2 
E 20 j 


(D22) 


?2 


_ *L e**R 2 
E 00 j 


The solution of equations (D21) is 


D /B 

j J 


1 

O m 
< >• 

1 

A A 

E /B 

' 

(_u W 
1 


(D23) 
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where 


? S = l (£ 22 d il - £ 12 d 20 )/(£ U £ 22 - £ 1?> 


AW 1 A p A /V A A A O 

Y j = 2 R j E (L ll d 20 ” L 12 d ll )/(L 11 L 22 " L 12 ) 


From equations (D24) and (D18) it is evident that 


(D24) 


-D 

Y 2 = 


-D -D 




Y 1 > y 4 = ‘ Y 3 


a E a E 

y 2 = - Y 1 ’ Y 4 


- Y- 


(D25) 


For j - 5 (r = r^ -= 0) equations (D12) give 


— — 


— ■„ — 



— c„ _ 0 0 


B c 


0 

22 


5 



0 ; 22 0 


£ 5 

= 

0 

0 0 0 


B 5 


0 

_ __ 


_ _ 


_ 


whence 



E,- = indeterminate 


(D26) 


(D27) 


Summing the four solutions of the form of equations (Dll) corresponding 
to r = r^ through r^ and making use of equatie»is (D23) and (D25), one 
obtains 


U 2 


11 11 



~D -'D ~D ~D 

V 2 


Y 1 - Y 1 Y 3 ' Y 3 



~E ~E ~E ~E 

r? 

o 

1 


Yl - Y X Y 3 - Y 3 _ 


B l e 


R. - 

1 c 


B 2 6 


-R - 
1 c 


V 


Ro - 

3 c 


B 4 e 


-St 


(D28) 
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Expressing the exponential functions in terms of hyperbolic functions, 
and taking into account the fact that U£ must be even in z , V 2 and 
vq odd in z , one converts equations (D28) to 


R. z R_z 

Bt cosh + B* cosh 

1 c 3 c 


R^z 


V 


B * sinh + Y 3 B 3 sinh — 


R,Z RqZ 

Y, B* sinh + y, B* sinh 

'1 1 c '3 3 c 


> (D29) 


where B* and B^ are new arbitrary constants arising from certain linear 
combinations of the previous ones. 

From equations (D27) and (Dll) the following solution corresponding 
to r = r^ = 0 is obtained: 

Ug ■ 0 


v„ * 0 


Since v o * E 5 = constant is even in z , rather than odd, the constant 
E 5 may be equated to zero. Thus the pertinent solution contributed by 
the root r = £5 = 0 is identically zero. 


Complementary solution not of exponential form . - The existence of 
a repeated root (see eq. (D19)) indicates that the homogeneous system 
has an additional solution that is not of the form of equations (Dll) . 
This additional solution can be obtained by inspection if one postulates 
that it has the following form consistent with the evenness of U 2 and 
the oddness of V2 and vq : 


u 2 = B 6 


v 2 =* Dgz/c 


v Q = E 6 z/c 


> (D30) 


J 
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Substituting this assumption into the system of differential equations 
(D8) with the ci2 u 0 term omitted shows that equations (D30) are indeed 
a solution of the homogeneous system provided that 


d 6 = 0 




where 


(D31) 


5 2 c*c 22 /d 20 


(D32) 


Therefore the additional solution of the homogeneous system is 


u 2 = B 6 


v, = 0 


v 0 = ?B 6 z/c 


> (D33) 


Complete solution of the differential equations (D8) . - Summing the 
particular solution (DIO) and the various solutions (D29) and (D33) of 
the homogeneous system, one obtains the following complete solution of the 
differential equations (D8) having the pertinent symmetry and antisymmetry 
properties : 


c 12 R l z R 3 z 

u„ u„ + B* cosh B* cosh + B, 

2 c 22 ® 1 c 3 c 6 


v„ = Y-i B * sinh + y~ B* sinh 

2 11 c 3 3 c 


-E R l z -E R o z 

v = Yi B* sinh + y B* sinh + £B,z/c 

U ±1 Cjj c o 


> (D34) 
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Evaluation of the arbitrary constants . - The three arbitrary constants 
B$, B^, and fig can be determined from the boundary conditions, equations 
(D9a) and (D9b) . These boundary conditions lead to the following equations 
defining B*, B^, and Bg : 


Q 11 

Q 12 

0 


* 

B 1 

— sinh 

u o 

Rjb 

c 

Q 21 

q 22 

Q 23 


B * 

— sinh 

u o 

V 

c 

Q 31 

q 32 

Q 33 


®6 /u 0 



where 


<11 


= R„ 


<12 


<21 


0 



0>35) 




Q 22 ~ Y 3 


Q 23 - 5 b/c 


> (D36) 


. .2 


R-b 


Q 3 1 - {d* x + 2R x [(e** + §“ \ e 22 )Y l + e 20 Y l ]}c ° th “T 


<32 


{d*i + 2i 3 [(e»* + 4 «$2>^ + *20 ^>“ th ¥ 

c 


<33 


d* + 2e** 
11 20 




2 sin6 

1 + 2 y 


(D37) 
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Relationship between F and u Q . - With u 2 (z) known in terms of 

Uq (the first of eqs. (D34)), equation (D7) yields the following relation- 
ship between the shearing force F and the relative shearing displacement 


where 



Gtb 

k 




* = 1 


1 f B ? *l bS | 

1 + 2 £ l u 0 C J 


- (*1 6lnh 
l “o 


R^b 


R 3 b U 0 


(D38) 


(D39) 


Equation (D38) can be used to compute the overall shearing stiffness of 
a single corrugation. 

A relative shearing stiffness ft , the ratio of the shearing stiffness 
(D38) to that of the same corrugation with continuous end attachment 
producing uniform middle-surface shear strain throughout the sheet, 
is given by 


ft = (1 + J 


(D40) 


Stresses . - The longitudinal normal stresses are identically zero at 
stations (tJ) , (l) , (4) and (E) . The non-vanishing longitudinal normal 
stresses v (2) along junction line (2) , obtained from the strains du 2 /dz , 
are given iir-dimensionless form by 



B* R..Z B* R.z 

— R sinh 1- — R_ sinh 

u 0 1 c u o 3 


(D41) 


The dimensionless middle-surface shear stresses, as obtained from 
table 2 and equations (D34) , are given by 
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T 12 C 

Gu_ 


u 2 ( z ^ 


- 1 


■\ 


- COS0 


-E B 1 - V -E B $ - V 

Y? — R, cosh - i - + y, — R, cosh — + C— 

1 u 0 1 c 3 u 0 3 C u 0 


t 23 c 

Gu. 


> (D42) 


u (z) ( v B* . 

- 2 £ 


f u„ 


-E "1 " R l z -E B 3 a R 3 z 

Y, — R n cosh -±- + Yo — R- cosh 
' 1 Uq 1 c '3 Uq 3 c 


u o 


"N / b* ^ 

4- a -f n A 1 C D — R 

R l z -D B 3 * 

pncVi -• — a/ ^ D /> a 

R,z^ 

aH A 

J 8 n6 ( Yl u 0 1 

COSfl T Y n J\ A cc 

c ^3 u 0 3 

c 

/ 


where the subscripts- 12 and 23 represent the plate elements 12 and 23 
respectively. 

From the rates of twist in table 3 and the displacement expressions 
(D34) the following expressions are obtained for the extreme-fiber shearing 
stresses due to the twisting of plate elements 12 and 23 respectively: 




where 



2 y cos0 W(z) 


w(z) = Yt -r- 


B* „ 
- R. 
U 0 1 


R 1 Z 


cosh + y 


"D 


B* . 

, - R, 
3 u 0 3 


— R n cosh 


r 3 z 




> (D43) 


(D44) 


The frame bending moments and associated extreme fiber bending stresses 
are zero at stations ( 5 ) , (l) , (4) and (5) . From the first of equations 
(E38) of reference 1, one obtains the following dimensionless expression for 
the extreme fiber transverse bending stresses or ' along junction line ( 5 ) : 


°©‘ 

Eu. 


1 - v 


c t_ 
2 k f 


1 + 2 j cos0 


1 + 


2 - 
* f 


v 2 (z) 


(D45) 
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(a) Cross section 



(b) Plan view 


Figure 1. - Configuration of trapezoidally corrugated plate 
considered in the present analysis. 
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(a) Point attachment at the ends of the trough lines. 




* 


/ 


(b) Point attachment at the ends of the trough lines and crest lines. 



//////A v- 



2/ 

'//A V////////7Z/ 


(c) Wide attachment at ends of trough lines only. 


(d) Idealization of (c) used in the analysis: Point attachments at 

the ends of the trough lines, and point attachments permitting 
longitudinal sliding at the junctions of the trough plate elements 
and the inclined plate elements. 


Figure 2. - Types of attachment considered at the ends of the corrugations. 




(a) Diagram of cross section. 


© 



Assumption regarding 
longitudinal displace- 
ments . 



(c) Component displacement 
modes for displacements 
in the plane of the 
cross section 


Figure 3. - Diagrammatic representation of assumptions regarding displacements. 




Figure 4. - Schematic representation of a 
single corrugation. 
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Figure 5. - Sign convention for t~, , t' and o’ . 

01 01 1 
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Figure 6 
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(c) Special case e = 0 


General and special cross-sectional geometries 
considered in the analyses. 
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(b/p)<t/p) 3/2 
(a) h/p = .1 


Figure 7. - Relative shear stiffness for the case of point 
attachments at the ends of trough line only. 



1.0 



(b/p)(t/p) 3 ^ 2 
(b) h/p * .2 


Figure 7 
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(c) h/p = .3 


Figure 7. - Continued 
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(d) h/p * .4 


Figure 7. - Continued 
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Figure 7. - Concluded, 
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Figure 8. - Relative shear stiffness for the case of point attachments 
at the ends of both the trough lines and the crest lines. 
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Figure 8. - Continued 
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Figure 8. - Continued 
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Figure 8. - Continued 





.00035 .001 .01 .1 1 2.5 


(b/p)(t/p) 3/2 
(e) h/p = .5 


Figure 8. - Concluded 
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Figure 9. - Relative shear stiffness for the case of wide 
attachments at the ends of the trough lines. 
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Figure 9. - Continued 
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Figure 9. - Continued 
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Figure 9. - Continued 
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Figure 9. - Concluded 
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Figure 10. - Concluded 
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Figure 11. - Special log-log plot of relative shear stiffness data for the case 
of point attachments at the ends of both the trough lines and the 
crest lines. 
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Figure 11. - Concluded 
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Figure 12. - Special log-log plot of relative shear stiffness data for the 
case of wide attachments at the ends of the trough lines. 
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Figure 12. - Concluded 
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Figure 13. - Comparison of shear stiffness for trough lines held straight (ref. 1) and 
trough lines permitted to curve (present analysis), for the case of point 
attachments at the ends of the trough lines only. 
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Figure 14. - Comparison of shear stiffness for trough lines held straight (ref. 1) and 
trough lines permitted to curve (present analysis), for the case of point 
attachments at the ends of both the trough lines and the crest lines. 
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Figure 15. - Comparison of shear stiffness for trough lines held straight (ref. 1) and trough 
lines permitted to curve (present analysis), for the case of wide attachments 
at the ends of the trough lines. 


























0 



Figure 18. - Variation of stresses and displacements along the length 
of the corrugation for a particular geometry with wide 
attachments at the ends of the trough lines. 
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Figure 19. - Dimensionless maximum-stress parameters for the case of 
point attachments at the ends of the trough lines only. 
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Figure 19. - Continued 





(a’lll) Extreme-fiber shear stress for h/p = .1 


Figure 19 . - Continued . 
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(b*I) Frame bending stress for h/p = .3 


Figure 19. - Continued 
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(b*II) Middle-surface shear stress for h/p = .3 


Figure 19. - Continued 
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Figure 19. - Continued 
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(c*I) Frame bending stress for h/p = .5 
Figure 19 . Continued . 
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(c*II) Middle-surface shear stress for h/p = .5 
Figure 19. - Continued. 
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(c*III) Extreme-fiber shear stress for h/p = .5 

Figure 19. - Concluded. 
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(a* I) Frame bending stress for h/p = .1 

Figure 20. - Dimensionless maximum-stress parameters for the case of point attachments 
at the ends of both the trough lines and the crest lines. 
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(a* II) Middle-surface shear stress for h/p = .1 


Figure 20. - Continued 



10 



(b/p) (t/p) 3/2 


(a*III) Extreme-fiber shear stress for h/p = .1 
Figure 20. - Continued. 
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Figure 20. - 


Continued 
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20. - Continued 







(b*III) Extreme-fiber shear stress for h/p = .3 
Figure 20. - Continued. 
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(c*II) Middle-surface shear stress for h/p * .5 
Figure 20. - Continued. 
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(a* I) Frame bending stress for h/p = .1 


Figure 21. 


Dimensionless maximum-stress parameters for the case of 
vide attachments at the ends of the trough lines. 
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(b*II) Middle-surface shear stress for h/p 
Figure 21. - Continued. 
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Figure 21. - Continued. 
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(c-I) Frame bending stress for h/p = .3 
Figure 21. - Continued. 
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(d*I) Frame bending stress for h/p 
Figure 21. - Continued. 
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Figure 21. - Continued. 
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Figure 21. - Concluded. 
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Figure 22. - Comparison of experimental data of reference 8 
(as cited in ref. 2), present theory, and the 
theory of reference 2. 
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